AD694099 


Distribution  of  this  document  is  unlimited 


AFTAC  Project  No.  VELA  T/7704 

ARPA  Order  No.  624 

ARPA  Program  Code  No.  7F10 


FINAL  ANALYSIS  REPORT 
ALEUTIAN  ISLANDS  EXPERIMENT 
OCEAN-BOTTOM  SEISMOGRAPHIC  EXPERIMENTS 


Prepared  by 

A.  Frank  Linville  Gary  D.  McNeely  Ronald  D.  Bauer 
R.  Fred  Howard  A.  Ronald  Mendina 


Terence.  W.  Harley,  Program  Manager 
Telephone:  1-214-238-3473 

TEXAS  INSTRUMENTS  INCORPORATED 
Science  Services  Division 
P.O.  Box  5621 
Dallas,  Texas  75222 


Effective  Date  of  Contract:  5  April  1967 
Contract  Expiration  Date:  31  October  1968 
Amount  of  Contract:  $1,054,991 


C 


ACKNOWLEDGMENT 


D  D 

K\ECnF™ 

This  research  was  supported  by  the  W  QQT  7  19S9  * 

ADVANCED  RESEARCH  PROJECTS  AGENCY W 

Nuclear  Test  Detection  Office  J  jjTM  t-nFriU TJ  Tb 

under  Project  VELA  UNIFORM 
and  accomplished  under  the  technical  direction  of  the 
AIR  FORCE  TECHNICAL  APPLICATIONS  CENTER 
Contract  No.  F33657-67-C- 1341 


31  July  1968 


This  dociimont  has  h^en  crpprovso 
for  public  reloase  and  sals;  its 
iflrt.iTbeflnn  is  unlimited  _ _ 


,0 


BEST 

AVAILABLE  COPY 


o 


Distribution  of  this  document  is  unlimited 


AFTAC  Project  No.  VELA  T/7704 

ARPA  Order  No.  624 

ARPA  Program  Code  No.  7F10 


FINAL  ANALYSIS  REPORT 
ALEUTIAN  ISLANDS  EXPERIMENT 
OCEAN -BOTTOM  SEISMOGRAPHIC  EXPERIMENTS 


Prepared  by 

A.  Frank  Linville  Gary  D.  McNeely  Ronald  D.  Bauer 
R.  Fred  Howard  A.  Ronald  Mendina 

Terence  W.  Harley,  Program  Manager 
Telephone:  1-214-238-3473 

TEXAS  INSTRUMENTS  INCORPORATED 
Science  Services  Division 
P.O.  Box  5621 
Dallas,  Texas  75222 

» 

Effective  Date  of  Contract:  5  April  1967 
Contract  Expiration  Date:  31  October  1968 
Amount  of  Contract:  $1,054,991 

ACKNOWLEDGMENT 

This  research  was  supported  by  the 
ADVANCED  RESEARCH  PROJECTS  AGENCY 
Nuclear  Test  Detection  Office 
under  Project  VELA  UNIFORM 
and  accomplished  under  the  technical  direction  of  the 
AIR  FORCE  TECHNICAL  APPLICATIONS  CENTER 
Contract  No.  F33657-67-C -  1341 


31  July  1968 


soisnes  ssrvloss  division 


ABSTRACT 


The  crustal  structure  across  Amchitka  was  determined  from 
two  in-line  reversed -refraction  profiles  located  northeast  and  southwest  of 
Amchitka.  Both  profiles  employed  10  Ocean- Bottom  Seismographs  (OBS) 
spaced  at  20-km  intervals  and  5-ton  shots  detonated  at  20-km  intervals  off  each 
end  of  the  OBS  line.  Traveltimes  were  corrected  to  a  reference  plane  4.  5  km 
below  sea  level  to  minimize  lateral  velocity  variations  associated  with  the 
large  changes  in  water  depth  in  the  area.  Depth  to  the  Moho  along  the  profile 
was  determined  from  first  arrivals  using  standard  refraction  analysis,  wave- 
front  methods,  the  time-term  method  and  model  perturbation  techniques. 
Calculations  show  the  Moho  refractor  to  be  dipping  towards  Amchitka  with  the 
root  of  the  structure  located  60  km  north  of  Amchitka  under  Semisopochnoi 
Island.  The  Moho  velocity  is  8.  1  km /sec  in  the  area. 

Nearly  all  of  the  first  arrivals  were  Moho  refractions;  therefore, 
digital  techniques  were  employed  to  enhance  the  capability  of  picking  secondary 
arrivals  ir.  order  to  obtain  upper  crustal  velocity  information.  In  general, 
secondary  arrivals  were  difficult  to  identify,  and  although  some  apparent 
velocity  measurements  were  obtained,  the  information  could  not  be  directly 
incorporated  into  the  crustal  model. 
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SECTION  I 


SUMMARY 


The  Ocean- Bottom  Seismograph  (OBS)  Aleutian  Islands  Exper¬ 
iment  was  conducted  between  5  April  1967  and  31  October  1968.  The  objective 
of  the  experiment  was  to  determine  the  crustal  structure  across  Amchitka 
Island  and  subsequently  to  determine  the  effect  of  that  structure  on  seismic 
traveltime  anomalies  observed  from  the  LONGSHOT  explosion,  which  was 
detonated  on  Amchitka  on  29  October  1965. 

Unit  and  equipment  preparation,  field  operations,  and  data 
analysis  were  the  tasks  performed  under  the  experiment.  These  tasks  have 
been  discussed  in  previous  reports: 

•  Operations  Report  —  a  detailed  description  of 
unit  and  equipment  preparation  and  field  operations 

2 

•  System  Response  Report  —  a  description  of  the 
determination  of  vertical  and  horizontal  amplitude 
responses  of  the  Ocean- Bottom  Seismograph  system 

•  Preliminary  Analysis  Report^ —  details  of  data 
quality,  basic  reduction  of  the  first-arrival  data, 
determination  of  charge  depths  using  hydroacoustic 
traveltime  recordings  and  bubble-pulse-period 
measurements,  compilation  of  preliminary  bulletin, 
and  a  model  of  the  crustal  structure  across  Amchitka 
obtained  using  standard  refraction  methods 

Final  results  of  the  data  analysis  program  are  presented  in 

this  report. 


Data  obtained  from  the  Aleutian  Islands  Experiment  were 
»  collected  in  three  phases.  Instrument  arrays  (OBS)  and  shot  arrays  (5-ton 

explosions)  were  deployed  along  a  line  N42°E  across  Amchitka  for  each  phase. 
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OBS  instrument  and  shot  locations  for  Phases  I  and  III  are  shown  in  Figures 
1-1  amd  1-2,  respectively.  These  two  arrays  were  designed  to  record 
ref  ction  arrivals  from  the  Moho  in  order  to  determine  the  crustal  structure 
alc.ig  the  profile. 

Phase  II  instruments  (Figure  1-3)  were  deployed  at  longer 
ranges  with  larger  spacings  than  Phases  I  and  III  to  record  seismic  waves 
penetrating  the  upper  mantle  from  the  CHASE-VI  explosion.  However,  CHASE 
VI  did  not  detonate,  precluding  direct  upper  mantle  observations  by  the  OBS. 

The  explosion  recordings  were  edited  onto  digital  tape  and 

played  back,  and  first  arrivals  were  picked.  Due  to  significant  changes  in 

water  depth  and  structural  environment  across  the  profile  (Figure  1-4),  the 

traveltime  data  were  corrected  to  a  reference  plane  (4.  5  km  below  sea  level) 

for  interpretive  purposes.  The  crustal  parameters  used  in  correcting  the 

3 

data  to  a  reference  plane  were  estimated  by  a  series  of  iterations.  Errors 
resulting  from  incorrect  velocity  assumptions  were  minimized  by  choosing 
the  reference  plane  to  approximate  the  water  bottom  on  either  side  of  the 
Aleutian  Ridge. 

Nearly  all  of  the  first  arrivals  obtained  were  Moho-refracted 
events;  therefore,  any  upper  crustal  velocity  information  must  necessarily  be 
derived  from  secondary  arrivals.  Digital  techniques  were  employed  to  en¬ 
hance  the  capability  of  picking  secondary  arrivals.  Traces  were  gathered  to 
produce  records  consisting  of  all  explosions  recorded  at  one  station  and  all 
stations  recording  one  explosion.  Static  corrections  were  applied  to  align  the 
traces  at  a  velocity  of  8  km/sec;  trace  amplitudes  were  equalized  to  reduce 
signal  variations  across  the  records;  and  deconvolution  filters  were  designed 
and  applied  to  remove  reverberation  effects.  Some  secondary  arrivals  were 
picked,  time-distance  plots  generated,  and  apparent  velocities  determined 
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from  least-square  line  fits.  In  general,  however,  secondary  arrivals  were 
difficult  to  identify,  even  after  processing.  In  addition,  a  positive  identifi¬ 
cation  of  the  same  secondary  arrival  from  opposite  directions  across  the 
profile  could  not  be  made.  Thus  the  velocity  information  derived,  while 
indicating  the  existence  of  intermediate  refractors,  is  qualitative  and  cannot 
be  directly  incorporated  into  the  crustal  model. 

Section  II  describes  the  four  methods  of  crustal  structure 
determination  used.  Only  first-arrival  information  was  used  for  each  method, 
and  the  average  crustal  velocity  was  taken  as  6.4  km/sec.  The  methods  are 
standard  refraction,  wavefront,  time-term,  and  model  perturbation.  Figure 
1-5  shows  a  comparison  of  the  resulting  model  obtained  using  each  method. 

The  structure  obtained  from  each  method  is  similar;  however,  differences 
arise  from  the  basic  assumptions  necessary  to  implement  each  method. 

The  standard  refraction  method  assumes  plane-dipping  constant- 
velocity  layers  across  the  profile.  Phase  I  and  Phase  III  OBS  instrument  arrays 
and  shot  arrays  were  designed  to  divide  the  structure  across  the  Aleutian 
Ridge  into  four  segments,  each  of  which  could  be  approximated  by  a  plane¬ 
dipping  layer.  A  reverse-time  mistie  showed  that  the  structure  under  Phase  III 

3 

instruments  could  not  be  adequately  represented  by  a  plane-dipping  layer. 

The  wavefront  method  (applied  to  the  structure  under  Phase  I 
instruments  only)  indicated  the  adequacy  of  the  method  there  and  showed  the 
need  for  denser  recording  along  the  profile. 

The  time-term  method  produced  distortion  in  the  model  associ¬ 
ated  with  the  root  of  the  structure.  Theoretical  calculations  verified  the 
distortion  which  was  attributed  to  a  deviation  from  the  time-term  assumption 
of  interface  planarity  within  the  cone  of  critical  refraction. 
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The  traveltime  model  was  determined  by  perturbation  techniques. 
Iterations  were  performed  by 

(1)  Selecting  a  beginning  model 

(2)  Calculating  theoretical  traveltimes  for 
the  model 

(3)  Modifying  the  model  based  on  differences 
between  calculated  traveltimes  and  observed 
traveltimes 

(4)  Repetitively  applying  steps  2  and  3  until  a 
best  fit  is  determined 

The  model- perturbation  method  produces  a  model  (not  neces¬ 
sarily  unique)  which  fits  the  data  by  considering  the  actual  ray  paths  involved 
and  is  not  restricted  by  the  various  assumptions  which  are  necessary  to  im¬ 
plement  the  other  methods.  In  addition,  the  method  utilizes  every  available 
shot/station  observation  along  the  profile.  For  these  reasons,  the  perturbation 
model  presented  in  Figure  1-5  (bottom)  is  considered  the  most  accurate 
representation  of  the  crustal  structure  across  Amchitka  which  can  be  obtained 
from  the  data  collected  in  this  experiment. 

Section  III  compares  traveltimes  observed  at  distance  ranges 
of  16°  to  24°  for  a  Jeffreys-Bullen  (JB)  crustal  model  and  the  Aleutian  crustal 
model.  Due  primarily  to  the  differences  in  average  crustal  velocity,  variations 
up  to  1.  5  sec  are  observed.  The  effect  of  a  dipping  structure  can  be  seen,  but 
effect  is  less  than  1  sec  for  the  model  used.  This  analysis  is  not  exhaustive 
but  does  suggest  that  a  significant  portion  of  the  traveltime  anomalies  ob¬ 
served  from  LONGSHOT  is  caused  in  the  upper  mantle. 
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Figure  1-2.  Phase  III  Station  and  Explosion  Locations 


Figure  1-3.  Phase  II  Station  and  Explosion  Locations 
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Figure  1-5.  Comparison  of  Models  Determined  by  Refraction,  Wavefront, 
Time-Term,  and  Perturbation  Methods 
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SECTION  II 
CRUSTAL  STUDY 


The  objective  of  the  study  presented  in  this  section  was  to 
determine  the  crustal  structure  across  the  Aleutian  Ridge  so  that  its  effect 
on  seismic  traveltime  anomalies  observed  at  teleseismic  distances  from  the 
LONGSHOT  explosion  on  Amchitka  could  be  defined.  Data  used  were  collected 
in  Phase  I  and  Phase  III  operations.  Phase  I  (Figure  1-1)  consisted  of  a  line 
of  10  OBS  instruments  south  of  Amchitka  spaced  20  km  apart,  with  10  ex¬ 
plosions  detonated  in-line  to  the  south  and  two  explosions  in-line  to  the  north 
of  Amchitka.  The  Phase  III  array  (Figure  1-2)  consisted  of  10  OBS  instruments 
located  north  of  Amchitka  along  the  Petrel  Bank,  with  eight  explosions  in-line 
to  the  north  and  three  in-line  to  the  south  of  Amchitka. 

Results  of  four  crustal  analysis  methods  using  data  from  Phase 
I  and  Phase  III  are  presented  in  this  section.  The  analysis  methods  are  re¬ 
fraction,  wavefront,  time-term,  and  model  perturbation.  The  refraction 

a 

analysis  was  presented  in  the  Preliminary  Analysis  Report  along  with  details 
of  data  quality  and  data  handling.  Those  results  are  briefly  reviewed  here  for 
comparison  with  results  from  the  other  methods. 

Because  of  the  shot  and  seismometer  geometry  used  in  the  ex¬ 
periment,  nearly  all  first  arrivals  were  Moho- refracted  events.  Consequently, 
a  study  of  secondary  arrivals  was  undertaken  in  order  to  obtain  upper  crustal 
velocity  information.  Trace  gathering,  amplitude  equalization,  and  deconvolu¬ 
tion  techniques  were  applied  to  the  digitized  data  to  aid  in  the  interpretation 
of  secondary  arrivals. 

A.  DATA  REDUCTION 

1.  Correction  to  Reference  Plane 

Since  all  explosions  and  receivers  were  not  located  on  the  same 
plane,  the  data  were  reduced  to  a  reference  plane  for  interpretive  purposes. 
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The  method,  described  in  the  Preliminary  Analysis  Report,  ^  reduced  the 
raw  traveltimes  to  traveltimes  below  a  reference  plane  by  stripping  off  the 
time  from  explosion  to  reference  plane  and  station  to  reference  plane. 

Figure  II- 1  illustrates  this  method.  Computer  programs  were  written  to 
process  the  arrival-time  data. 

The  following  parameters  are  input  to  the  correction  program: 

Vj  —  water  velocity 
V 2a  —  subwater  velocity  (low) 

V23  —  subwater  velocity  (high) 

Vn  —  refractor  velocity 
RP  —  reference-plane  depth 
DP  —  decision-plane  depth 

VRed  —  velocity  used  to  compute  reduced  traveltime 
t  ,  1 .  e  •  , 


where 

t  =  traveltime  corrected  to  RP 

A  =  distance 

The  output  of  this  program  is  then  used  to  generate  CalComp 
plots  of  the  raw,  corrected,  and  reduced  traveltimes  vs  distance.  In  addition, 
a  third  program  computes  least- squares  line  fits  to  both  the  raw  times  and 
corrected  times. 

2.  First  Arrivals  and  Secondaries 

The  refraction  model  developed  during  the  preliminary  analysis 
was  based  on  first  arrivals.  Those  arrivals  were  grouped  into  three  categories 
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Pn,  Pg,  and  shallow  Pg  arrivals.  Based  on  first  arrivals,  portions  of  only 
two  refra  tors  above  the  Moho  could  be  mapped.  Therefore,  additional  upper 
crustal  velocity  information  must  come  from  secondary  arrivals. 

Attempts  to  pick  secondary  arrivals  from  the  field  records  did 
not  produce  meaningful  results;  consequently,  additional  processing  was  per¬ 
formed  to  enhance  the  capability  of  picking  secondary  arrivals.  Digitized 
explosion  recordings  (Figure  II-2)  from  Phases  I  and  ITT  and  the  northern  part 
of  Phase  II  were  used  as  input  data  for  TIAC  processing,  which  comprised 

•  Trace  gathering 

•  Trace  amplitude  equalization 

•  Static  correction 

•  Paper  playbacks  (filtered  and  unfiltered) 

The  initial  trace- gathering  sequence  consisted  of  gathering  the 
pressure  and  vertical  traces  for  all  explosions  recorded  by  one  station  to  form 
one  output  record.  The  two  horizontal  traces  were  gathered  in  the  same  manner 
to  form  a  second  record.  Trace  amplitudes  were  equalized  and  static  cor¬ 
rections  applied  to  align  the  traces  at  a  velocity  of  8  km/sec  (Figure  II-3). 
Equalization  was  performed  by  scaling  the  traces  so  that  their  mean  absolute 
values  over  the  equalization  gate  (9.6  sec)  were  equal.  This  procedure  gives 
record-to-record  equalization,  as  well  as  trace-to-trace  equalization.  Using 
the  records  of  the  first  gathering  sequence,  the  traces  were  gathered  to  form 
records  consisting  of  ail  stations  recording  a  single  explosion  (Figure  II-4). 

A  second  set  of  records  was  made  by  applying  static  corrections 
to  individual  traces  in  order  to  vertically  align  the  first  motion  of  all  traces 
(Figures  II - 5  and  II-6).  After  each  gathering  sequence,  paper  playbacks  of 
all  records  were  made  and  frequency  filters  applied  to  further  enhance  the 
signal  quality.  Examples  of  the  secondary  arrivals  are  shown  in  Figures  II- 7 
through  II-  17 . 
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Figures  II - 7  through  II- 10  are  records  showing  all  explosions 
recorded  by  a  station  (pressure  and  vertical  traces).  The  predominant  arrival 
after  the  initial  P-phase  is  the  arrival  corresponding  to  one  air/water  inter¬ 
face  surface  reflection  in  the  water  layer  at  the  receiver.  Only  in  the  case  of 
all  explosions  recorded  by  S10  (Figure  II- 7 ,  top)  could  an  arrival  be  followed 
across  all  traces.  Most  of  the  later  arrivals  appear  to  be  emergent. 

Figures  II- 11  and  11-12  show  examples  of  the  horizontal  traces 
for  all  explosions  recorded  at  a  single  station.  Stations  S10  and  S9  show  the 
best  arrivals,  but  only  two  or  three  picks  could  be  made.  Records  of  S2  and 
SI  show  the  cyclic  nature  of  many  of  the  arrivals.  Also,  it  can  be  seen  that 
the  moveout  pattern  is  not  consistent. 

Shown  in  Figures  11-13  through  11-15  are  examples  of  all  stations 
recording  a  single  explosion  (pressure  and  vertical  traces).  As  in  the  other 
examples,  the  secondary  arrivals  are  emergent,  with  the  number  of  picks 
scattered  from  record  to  record.  Figures  11-13  (bottom)  and  11-14  (top)  show 
the  air/water  interface  surface  reflection  at  the  station. 

Phase  III  recordings  are  illustrated  in  Figures  11-16  and  11-17. 
Figure  11-16  shows  the  pressure  and  vertical  traces  for  S31  and  S23,  and 
Figure  11-17  shows  the  horizontal  traces  for  the  same  stations.  All  secondary 
arrivals  for  Phase  III  were  picked  from  recordings  showing  all  explosions 
recorded  at  one  station.  Records  showing  all  stations  recording  one  explosion 
were  of  little  value  because  of  high  noise  levels  or  instrument  problems  at 
many  stations.  Secondary  arrivals  from  Phase  II  were  insufficient  to  produce 
meaningful  data. 

After  all  secondary  arrivals  had  been  picked,  the  times  were 
converted  to  GCT,  punched  on  cards,  and  processed  through  the  correction 
program  using  the  following  parameters: 
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V ^  =  1.5  km/sec 

V2  =  2.1  km/ sec 
^2B  =  5.  0  km/ sec 
Vn  =  7.0  km/ sec 
RP  =  -4.  5  km 
DP  =  3.  5241  km 
VRed  =  8.  0  km/sec 

Output  data  from  the  reference-plane  program  were  divided 
into  four  groups  for  plotting: 

•  Pressure  and  vertical  arrivals  —  all  explosions 
recorded  at  each  station  (Phase  I) 

•  Pressure  and  vertical  arrivals  —  all  stations 
recording  each  explosion  (Phase  I) 

•  Horizontal  component  arrivals  (Phase  I) 

•  All  secondary  arrivals  (Phase  III) 

Plots  of  the  raw  times  and  corrected  times  (Figures  11-18 
through  11-25)  were  made.  Comparison  of  the  two  sets  of  plots  shows  that 
the  largest  effect  of  the  reference-plane  correction  was  to  lower  the  points 
by  2  to  4  sec  (which  is  essentially  the  traveltime  through  the  water  layer  from 
shot  to  the  water  bottom). 

Reduced  traveitime  plots  of  the  four  sets  of  data  are  shown 
in  Figures  11-26  through  11-29.  Lines  connecting  the  points  which  represent 
secondary  picks  across  the  records  are  shown.  It  is  apparent  that  some  picks 
\  represent  multiple  reflections  in  the  water  layer.  Water  depth  at  station  S2 

is  1.44  km,  which  gives  a  2-way  traveltime  through  the  water  layer  of  about 
2  sec.  The  three  sets  of  secondary  arrivals  picked  from  the  horizontal  in- 

f 

struments  at  S2  (Figure  11-28)  show  a  time  separation  of  abour  2  sec,  indicating 
that  the  two  latest  sets  were  multiple  reflections  in  the  water  layer.  Similarly, 
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the  two  latest  sets  of  secondary  arrivals  (Figure  11-29)  picked  from  station 
S31  (water  depth  of  1.58  km)  show  a  time  separation  corresponding  to  mul¬ 
tiple  reflections  in  the  water  layer.  Water  depths  for  stations  S5  through  S10 
(and  the  southern  shots  of  Phase  I)  indicate  a  2-way  traveltime  of  6  to  9  sec; 
therefore,  the  secondary  arrivals  occurring  in  the  6-  to  9- sec  interval  follow¬ 
ing  the  initial  arrival  cannot  be  attributed  to  multiple  reflections  in  the  water 
layer. 

Least-square-error  line  fits  (Table  II- 1)  were  made  for  each 
set  of  arrivals  having  at  least  four  picks.  Line  fits  for  Phase  I  explosions 
recorded  by  individual  stations  (pressure  and  vertical  traces)  gave  apparent 
velocities  ranging  from  6.89  km/sec  to  7.96  km/sec.  Velocities  from  pressure 
and  vertical-trace  picks  of  gathered  records  (stations  recording  explosions) 
ranged  from  5.70  km/sec  to  6.65  km/sec.  Individual  line  fits  for  secondary 
arrivals  picked  from  the  horizontal  traces  gave  velocities  ranging  from 
5.  13  km/sec  to  7.53  km/sec.  Due  to  the  data  quality  of  Phase  III,  only  eight 
line  fits  for  exp  osions  recorded  at  stations  could  be  determined.  The  velocities 
ranged  from  6.61  k  n/sec  to  7.43  km/sec. 

In  addition  to  the  individual  line  fits,  the  secondary  arrivals 
were  grouped  to  determine  average  velocities  (Table  II-2).  Phase  I  arrivals 
from  all  explosions  recorded  at  one  station  gave  an  average  velocity  of  7.42 
km/sec.  Phase  I  secondary  arrivals  for  all  stations  recording  one  explosion 
gave  an  average  velocity  of  6.23  km/sec.  The  average  velocity  for  Phase  I 
horizontal  arrivals  was  6.67  km/sec.  Phase  III  arrivals  gave  an  average 
velocity  of  7.  04  km/sec.  The  secondary  arrivals  indicate  the  existence  of 
intermediate  refractors.  A  single  intermediate  refractor  with  a  true  velocity 
of  about  6.  8  km/sec  and  dipping  toward  Amchitka  could  explain  most  of  the 
measured  apparent  velocities  under  Phase  I  shots  and  stations.  It  was  im¬ 
possible  to  positively  identify  the  same  secondary  arrival  from  two  directions 
across  the  profile;  this  is  necessary  to  determine  true  velocity,  dip,  and 
depth  of  the  refractor. 


II-6 


science  services  division 


Results  of  this  study  show  the  difficulty  of  identifying  secondary 
refraction  arrivals,  picking  their  onsets,  and  identifying  the  same  arrival 
recorded  from  opposite  directions.  The  techniques  employed  here  —  trace 
gathering,  static  correction,  and  trace  amplitude  equalization  —  optimize  the 
recorded  data  for  interpretation;  however,  the  results  indicate  the  importance 
of  utilizing  shooting  and  recording  techniques  which  allow  recording  first 
arrivals  from  each  interface.^ 


Table  II- 1 

LEAST-SQUARES  LINE  FITS 

Phase  I,  Shots  Recorded  at  Stations 

Velocity  Intercept 


Station 

Explosion 

(km/sec) 

(sec) 

Slope 

SI 

E34,  E3,  E5,  E6, 
E30,  E31,  E32 

7.32 

2.  54 

0. 1365 

S2 

E5,  E6,  E30, 

E31,  E32 

7.68 

3. 07 

0. 1302 

S5 

E2,  E3,  E4,  E6, 
E30,  E31 

7.27 

3.63 

0.  1373 

S6 

E2,  E3,  E4,  E5, 

E32 

7. 57 

3.57 

0.  1321 

S6 

E6,  E30,  E31, 

E32 

7. 96 

3.31 

0. 1256 

S7 

E34,  E2,  E3,  E4, 
E5,  E6 

7.24 

4.40 

0. 1382 

S9 

E34,  E3,  E4,  E5 

6.89 

1.70 

0.  1452 

S9 

E3,  E4,  E5,  E6, 

E30 

7.42 

2.75 

0. 1348 

S10 

E34,  E2,  E3,  E4, 
E5,  E6,  E30,  E31, 
E32,  E33 

7.46 

3.41 

0. 1341 

S10 

E3,  E4,  E5,  E32, 
E33 

7.41 

6.68 

0. 1349 
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Table  II- 1  (contd) 


Phase  I,  Stations  Recording  Shots 

Velocity  Intercept 


Station 

Explosion 

(km/ sec) 

(sec) 

Slope 

S2,  S5, 

S7 , 

S8 

E2 

6.64 

2.  39 

0.  1505 

S2,  S5, 

S6, 

S7 , 

S8 

E4 

6.44 

0.94 

0.  1552 

S2,  S5, 

S6, 

S7 , 

S9, 

S10 

E4 

6.  24 

8.89 

0. 1603 

S6,  S7 , 

S8, 

S9, 

S10 

E4 

6.47 

2.58 

0. 1646 

S2,  S5, 

S6, 

S7 , 

S9, 

S10 

E5 

5.69 

2.42 

0.  1758 

SI,  S2, 
S10 

S5, 

S6, 

S7 , 

S8, 

E5 

6.  15 

11.84 

0.  1627 

SI,  S2, 

S5, 

S6, 

S9 

E5 

5.70 

12.87 

0.  1754 

SI,  S2, 

S5, 

S6, 

S7 , 

S8 

E31 

6.65 

1. 58 

0.  1504 

S7,  S8, 

S9, 

S10 

E31 

6.09 

-2.51 

0.  1642 

Phase  I,  Shots  Recorded  at  Stations  (Horizontal) 


SI 

E3, 

E4, 

E5, 

6.  38 

1.21 

0. 1568 

E6 

S2 

E2, 

E3, 

E4, 

E5, 

E6, 

E30, 

E31 

7.44 

4.50 

0. 1344 

S2 

E2, 

E3, 

E4, 

E5, 

E6, 

E30, 

E31 

7.53 

7.48 

0. 1328 

S2 

E2, 

E3, 

E4, 

E5, 

E6 

7.21 

8.69 

0.  1388 

S9 

E5, 

E6, 

E31, 

E32 

6.61 

3.40 

0. 1513 

S9 

E3 , 

E4, 

E5, 

E6 

6.36 

10.86 

0. 1572 

S10 

E3, 

E4, 

E5, 

E6 

5.  13 

6.68 

0. 1948 
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Table  II-  1  (contd) 

Phase  III,  Shots  Recorded  at  Stations 

Velocity  Intercept 


Station 

Explosion 

(km/sec) 

(sec) 

Slope 

S31 

E16, 

E22 

E17, 

E18, 

E19, 

6.  95 

12.  14 

0.  1438 

S3 1 

E16, 

E21, 

E17, 

E22 

E18, 

E20, 

7.  18 

14.  58 

0.  1392 

S30 

E16, 

E 17 , 

E19, 

E20 

7.43 

13.95 

0.  1345 

S28 

E15, 

E16, 

E17, 

E18 

6.61 

21.84 

0.  1513 

S23 

E17, 

E21 

E18, 

E19, 

E20, 

6.79 

3.61 

0. 1473 

S3 1 

E15, 

E17 , 

E19, 

E21 

6.96 

16.  10 

0. 1436 

S28 

E15, 

E20 

E16, 

E17, 

E19, 

7.  33 

7.76 

0. 1364 

Table  II-2 

AVERAGE  VELOCITIES  OF  SECONDARY  ARRIVALS 

Velocity 
(km/ sec) 


Phase  I,  Shots  Recorded  at  Stations  7.42 

Phase  I,  Stations  Recording  Shots  6.23 

Phase  I,  Shots  Recorded  at  Stations 
(Horizontal)  6.67 

Phase  III,  Shots  Recorded  at  Stations  7.  04 
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Figure  II- 1.  Method  of  Reference- Plane  Correction 
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Figure  II-2.  Example  of  Digitized  Field  Data  from  OBS 
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Figure  II-3.  Gathered  Traces  of  All  Shots  Recorded  at  One  Station, 
Aligned  at  8  km/sec  with  Trace  Amplitudes  Equalized 


II- 12 


science  ssrvlcss  division 


E5 


STATION 

NUMBER 


2 

3 

CO 

</> 

2 


5 

6 
2 


•10  SEC- 


2 

5 

to 

2 

CL 


E15 


& 

2 


S  31 
30 
28 


27*++ 


26  m 

25 


S  23 

"sir 


pywwy^  ywyvs/v  *A/y<*<*  w*i4<«rfv  v*A*y*/W  «vvwrv*A*f*/wvw/v< 
s/^/AA  yv*r|yvW<  ,aam4+vu  *ys>Wv^  ntyu|>vM* 


:) 


28 

27 


25 
S  23  i 


Mh*n 


mm 


■o- 


'mw* 


>  ^  ^  S "/  Y  >  *,  V  < .  t,  y|\  S .  >  ti  >,  %  \  f  n  *jwf- 


hW 


U^/K+441 


m 


26 


d 


)w4w 


<  >4*M* "  4^  ■  *Wl)M  nM*  4*H‘  >M4mI 


r  A/na  a/w\m 


*s(M  ’SMi/W  ^ovw  w/V^vaAAsAw^^ 

xlA,  '  ^Ax^AvvA|v\AAA,4v\/KW\^AAoAA^ 


■10  SEC- 


Figure  II-4.  Gathered  Traces  o£  One  Shot  Recorded  at  All  Stations, 
Aligned  at  8  km/sec  with  Trace  Amplitudes  Equalized 
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Figure  II - 5 .  Gathered  Traces  of  All  Shots  Recorded  at  One  Station, 
Aligned  Vertically  with  Trace  Amplitudes  Equalized 
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Figure  II-6.  Gathered  Traces  of  One  Shot  Recorded  at  All  Stations, 
Aligned  Vertically  with  Trace  Amplitudes  Equalized 
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&  *  ^  v  ’  X*  v  f  .  ^  ;  v  »  ^  *  * 

3  - — n^wn a^v'WWVM'v  _  jV_  W  /;,  lVa  y\\  A  ■  /  ’  ; 

4  - . — yy  /  a’/  / /•’* <v  *  V  >' 1 

|  5  >_ - , \yV A 7  >a  ,J  •  V/  ;' 

S  6 - ^w'..  ,v~*v/*v * AWvv('vVy'(^vA^;  /v\  _VA  V  .y^AX  Wy\A'\ 

ss  30 - v-  v^vyvv.v.'w' .  v.v.v-  /.vy-y  /X  ■•  %  y 

31  - \  ;^vw/waw;  :»\v/y;'  '.X^vy  y 

32  - - AvAV\/'/sa/^,V,V'  /A'/.’Ay 
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Figure  II-7.  Gathered  Traces  of  All  Shots  Recorded  at  S10  (top) 
and  S9  (bottom) 
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E  34  -  •  J  ,/\Ay, \.  ^/Aa/V/A^A  \  V/  '■ 

2  - 'AAflvv — — w  •  *  wjwmVlvw 

3  - - wW  .,v  ,  (  I..  fj  ^V\A,  ii  ,  ;#A  *  .  f  , .  J 

, — svvww; — w^tj' •:*[•  j.i  ■  ■■  ■ 

6-W/*"  /r.ftt  J1 ....  >.M  ')  /,.  V.;L,  *  -J  '  in'v 

Y*34p- .  ■I-.IWWi/WiMA/AW  ■  ■l-xA'IWW  /7  j(  I ,  ',\\ 

2  . . !  ^VV^/VWVsWVVv^>^s»<^s^»i\/WA/^Vj/\/WWw|s/^\^X< 

3  — - — k^/V^Y/ww^avvw^  .  y»A>^y\/<Ww\Wvwyvv'/A'\/\wvwsvw\AWWvVv 

4  . .  — —/ \,  ’  /^V/A^^w^^^-^\Vv,y'-^w^A^wv^w^^A/^V>>^vVAAA//^Yv^A^wr 

’  \"  .'  ',  'y\V  '  V  .  y/  1  .  ‘I  •  v  ' 

30  '  1  •  '-^. ‘.V.V<V/^-^. l-^V.>wA)  ’ — >w\ ywv  V* — V; 


\Z^W/< 


(W  ,  _ V/'^^WA^, \". Vv* j\ V47a\vwa^ - ■wv\‘/* 

\A<vAAw/VNrAy^A\'\AA/'.'Ww>AA^v\  /WV.VwAW//,/V'-~-"v\v^~>V\«' 
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(,  ^/Vk4>/VN'vk^^WVl 


.''.V<v^ - - - V/Wv'  f  *:*•  *,  •  <vVwV  •“‘jVs  ^  */{fA  V.,.v\\v‘  aV“  •/ 

^'^Vvw - vw,  •  ^  v  a,AVvxWv,^/\'<V'v‘^V^vV  y/vj  ^ 

v^Wveww www^/  ,VA  #  ^WVV  A^A/V\  fvV :  /,  .i^vwM  A  v\  ;  A 

v,  ^ l^V,  _  V\  \\  4" ! . 

\  ’  “vv/Vn-aK  ,•  w  (  'v^/^vW“ \z"'\Av\,*yV/\  ;'Vv*v‘vVw'~  ‘/  VsA 

^A’M^>4^wvwwvA , I  \^i‘  Vv’  VyWV;V 

vSVvWwt-\M\'  4  A  Wav^vV/^V  r>:\  yvVA  V i  ■•/ ,v: 
-Vw^wVAN^^^A^YV^JV\^wA''^/y^  v  wjIvwAAVv  v/vA*j 

^/WVM  Aa^\a^a>n>«A^AW/  ,  1  ^  1 1\"  iWaAMmA  A  a/  nA  1  V  r  I'T/wAh  iaA^/v/ 


31  - - - — . 

32  - |WWV^MAa^  i\V\  .  i 

34  I” - -vsMyvvvv^  — ^A'A^vAW/Js^ 

3  1 — - — ■'vv'-aV  Jaa^7\A/VSi/'aA/  a*/v\,V  '■1  />  V  V\ ,  . W / „/k  \ j, lv\  ttyv\A 

4  - ,/WWV^ — A/vw\V^(  N\f4'v^AA^VWW\\\',^A/l'  lyW  \/  ^  V\'  \\/  v 

5  i  ■  y^v^^VVA\VwwW>vA//  _  jV/A\ww\^^\^\V///.-vVVyvWA\Vk\\w“'.''‘//' 

4  . —  ',W^Wv\ i/A/ws^-vCV-^ArWV*  \'  _ \\ \V^k Vi\vi\\^>''^Vv^v«^'i‘  '/v — v\  '///■' A’j'vwv 

30  ■— —  -  *  ■  ■  '^‘ywv-vVyv^'/.’^'vAV  7/'v  V.'.Vv, 

31  ■—*  ■  v'^VA/V,/WVv-w\rw^t‘7<\\V^‘  VW  \i4\\,»vv\(  \\\\ 

32  — - - —^AAA^W\^V\,Vw^VW\AfN^A^vv/yvwvry//(\^vw^>y<>kvc’A'^WVv\7.'.;.''A'w~v,.\v^V.'  '■'.' 


4\  /v  f,n 

*  aMVA 


1  i 


/;r,v^ — \\\ 
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Figure  II-8.  Gathered  Traces  of  All  Shots  Recorded  at  S8  (top)  and 
S7  (bottom) 
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^•vV^sAaw>^1'aV'  \  V\’ 

'  \/WV'\^-vV IJ '  \J\A/~*/\/  \<\A\  ’  ^wA  \  A/V*/ ;  V* 
^  ^^\Z^vVs,Vv  VAAA\  tVv*V^1'  t  * v 

A^>r\A/VvWVv^,v»' — VV\ \  V ,  a^'AVVv — W1 

A~\yv* - V\/\f '  ^'Wv^«aAV^  VA  VaA  1  ^ '  i, 

A,,^v»\/v/</V~ — ^A/vyw1  ^/v"  ^V'mA/a1,^  i/vv\A^/v 

a  Vv^'1\'V^V^ 

A^W\^^aAAaVV^^w Vv v  w  ¥  V\  *  yvV  '  /'vV 

VW^\M^VVvWv  «>/ AvV^S  ,^\Vv'i'1'/1'^v^/  /W*  ^V 
<\ — vWVWVWW^ \fX~J\sJ\ \ " w * a V'A'A '  V\A f  \A 

^/rW\Z<l\M'V^w,V/A‘;'M'.V\v->,  \/\_ '  .Wa'.'sa^ 


’  v  ^  V 

w  v  *  \v 


/V\-Wv^vVvW^k  .VWA  A^V.  y\\  VA\',“">\'\"'-/'>-V 


\V\7  V7  V.  ,W/.*-v’  ~'_  .  v'^V.  \'^ 
aaw^\\v^'v\',>V^;,',' 


-10  SEC- 


2j^— ' n-^-a  ^'.VNAA/WWVvv/ A^y/VyAA*7,*V  uVm^'^v) 

- - - / y-  /aaAwiVw^  V  >.'/~^a 

4  h - — \.  Ai  ,/v-V, V'-Vvrv^J\/'w V'w\"v'a  - 

tt  5 * - — ♦■*  ’  //a~Wvw^--~7  'i t\\rt VAWwi'://iM ’ ( \\7a\ V-n/^ 

gj  jo  -'■'^ — tW'/^'/v> '  f a/'*'>-4aa\  v//\A/// /.W,''^’/,’/ v/.^vy^AV  _  W 

31  '-^’v'^'U/ ,  \'V^W^A^WVv  V  /\Av  - Vv^ ^  W/^ f  *  A  V  VW\M/V\, 

32  '>^-s^^\V'/AA^w^VsyW\A^  'Mw^ 

*  k  E  33  |WVV^^W\np^,  Vv\  Vv'AyAA'^'^yA^V 

S  ^  E  34  ---*■■  4-^WV — v^a^a/v^V/i^v^A^/a*}^// /j'/%^'yvi>v/'jv\^V\^v'i, 

2  — - >^\nW/vvv+*sW^^ 

3  - -^V\A ?/  Vvvy^^^.v-^  rj  .**  >a/\V  *A‘,"  \  <*.  A\  '^V 

4  ^ wv  ,  /yv\^  'J/y  /*,  }s*Sy\rX 

o  5  -■-■  ■-"/'  • '  /.v/aaV 

j*  6  '— - *A ^//vV^vW/vWw'A^W^\'A^^\A^v^V/\>'X!>V(V'jy,  ^ 'M 

>  30  ‘  — — — —VvV\' ( V« ‘ /vVs^w'v v " A/ v. ‘ vW\A ' . V^-v'  /  '  ’,'S,\"  ‘/-* /J 

31  y^J^\r^Wvv\  V  AW,  \‘»/'/-\A'Vv 

32  \,^/^',\/sA,v/A^V',AA,,"^A^ A7*/\a/^>A/v  'v  /*AAvlW^ 
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Figure  II-9.  Gathered  Traces  of  All  Shots  Recorded  at  S6  (top)  and 
S5  (bottom) 
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4  |"i  »>i i  ^  ^*v  1 V~*v*Aa^*va/NA/s/s/\a/>/\/\AAaA'*"*%aA/V\//^/V^^ 

5  \W'/|^/V^WV^^^\MWW^,'/SA^^'vVvA,/>^vv^/' 

6  h- — >»w/iwv  v  \A^^|^^vv^^-^✓vv^/wv^y'l  Vy\, "  v  V/^v^a'a  V  vv 

g  -^v'/v/w\A/Yt/NV^^^^AA^y*/VNA^^vl‘w,^V~^VVV,'v\'V\^/VVv\/v<\1A, 

1  ■— -^'^nAVn/V^VV^^'^AJVV^^'  V'WV/^fV^v\'/°\A%V^,/V 


3  . - -V^-A/AV^  Vv'W'\*/-'AA/\/\  A’j  \*'/vV\(\, \ 

4  >  -  ■—  —  -  — "*^vt  V— \‘  <^✓«^V'/,,^vs/vv*/\/'v^ VIA^'/*v'/W\aA^-iW\/\/\/*'<>/WVvA 

j  u — - — V^/*  r/-^/A^^AV^^/A'-l*VirsA/'-^,'^AW /vV'Vvn^-''- 

4  P-^-^^AV^^^VvVwVv's~A~^^/V^\/^W^V^v^VA^.\Awvr y/'A/sXV 

31  <* — AA^^,y^WNW^vv  *v/  Vaw*W  ✓vV^V\ 
-  *  _  -  r  ^  _  V  »  .  .  / 


E  33 


V\-A vr^*W/ 


2  K-^^v/V'A/A/v\/vA'sA/VVv^'“''V*'^'W'-^,‘,'-,v\/\/W^/WV/^/''''rsA’N'''/^'v\/^' 


^  j  »-■— ■  «>,  .i  ^V^A/\Avw»^W\AA/\)l^A/VVylv>,^^V\;VWyvA/><s^y~wv 

&  6  wv\/vvy“v>^N''A/V/s*»\A'v*^/vA'  A^/\\/\/''\/^VV 

^  ’O  wv — 'N 


E  33  ’'Ma-vA/^v*>, 
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Figure  II- 10.  Gathered  Traces  of  All  Shots  Recorded  at  S2  (top) 
and  SI  (bottom) 
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SHOT 

NUMBER 


3  ~ 

4  ~ 

5  " 

6  - 
30” 


vVv\A Vvv^VvA  'WWW'WWW' 
*  \\ ’vVW,  \  \  Av VV/yW ^Vvyvv'-^,^V  Wwv-v- 

//.  /n.\’  '/V//  V 

^Avi^Ar^ /  A  W \AVw\V  '  V( •^-vf //a 
Avv^~rA^'jA/^A'w  \  \ ' 


JJ  ■  I  . VAV>^/AWVVWV^^A^;VV,AVV'JWVVf  JW^V\|  yw  J  /v 

32  ■-*«-—.-  -  v^Av^^^^A<\^v^AVS,Vw\V^/jVNAAA/\(\>y\AW*A,\//W  A'lW/ 

£  ■  I  ■  i  >  ■  Av'v1 WWN^A  l\/V^^»i  1  1  ^V»*  ■  —  ■—  ■-■■*^»/ ■■  '  ^  1  1  ■  .-<.•  -—  I  - 

3  i  i  -  ^■~AA/^AA^/V\/WAA/^t,lA/^^WVvw^W>v^'WV^,,vsvs^^^ 

4  >,  -  ■■■  ■✓vA,r'AA^^w^‘N^*/,^AVA\A\;\A>-^vVs/ A* /v"iM'Av<'v/V^ 

5  — i — ■ - w,Va\Vw\V — v\/'vv\y\  ‘  .-vaVv' .  *  _  VW/,  Y 

$  . -  *  AW;  VvVWv/vwv  \  ‘  y\  V .  Vv\  V**^  \  \ \  V\,  VW  AW  W‘-  ' . " , 

30  —  -  - - ^vYjV^aa/m'.V/.^.'^^w’/.vaaVwVwv 

31  - ;\\v\VA/\V 

32  .-  .-  —  ~^AW*<vvW/^^Vi^^tW><A^vv^.\VtfA<'l\‘  _  -s V/a 

£  33  - - - ^^\vwvv\^Vnv'>^wV^/a'>\\  w',\Aw^Ov\W\’  V  \-v 
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E  34  — — — — w^v^/wwwv'1  yA-^ ,  V/A,\ 1 

3  — — - - ■ — 

4  - — ^'wv-A^-«v;AA\Ay  WwV  /\V.'a^a\v 

Hj  5  - - - vW/^MW^AAV/k \/W?‘  k\V  I  V  "  Vj’lA^.*V  *  '  '  N 

6  — — V.  VW*v  VwWwW\,  VNAA/V  \  WvA  VAaW  V  ^  **  */ 

30  - ^;,v,'(aw/aa^w;a];  ///  'Uv^’jVvM'A'y;; 

31  - AWl^VWVv^W^W,  V^W/WWtVVi  „  AM  -V  • 

32  A^V.^',VvV/Vwa/A^  *  v  ) 
cq  E  33  |r^AAAWlMAV^ 

E  34  ^W/V  fWVKv — \W v'->/V^p^W\jr^^ 

3  - ^WvW/i\VW/  / 

4- - ■/Wv>~  -  v«..  v  wwyw*WP\  ,  'SA  i/l'"  *  A  MV.V .  V,  V  ’  M* 

5  — V'^.VVvW^ W/WVa'^a^/  "  _’  ‘.1 7  * '  V,  V  ’  r  ’  ;l  s  _V  ■ 

6  - /,\WM Aav^vAaA/'vW/.^^M  w’  '  \r\tjJ'j*JV]  * '‘  '  X’  .V 

30  - v7.WvV//A//,\\\W vaHWAV  'V  V/  \\W  J\'^\'}^'\\ 

31  - - - A*A'AAAAAAvVvW//vrvVW/A^Xs/y^AAV  V'V  V.V 

32  — - ■'AMa/sA* - vWWW/,y_Vv Wvyr//\V,,/\AA^/,\‘,AA/'/  ^  '  *< 

E  33  v;\'  ‘  7  V,y^/A'  ’/a Vi’;  j 
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Figure  II- 11.  Gathered  Traces  of  All  Shots  Recorded  at  S10  (top)  and 
S9  (bottom),  Horizontal  Instruments 
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Ww\ "w\ ^ V\v^s/\vv  <— 'ATA* 


6  VWVWVWWVW 


5  MAV, 


32  MAAMW, 


Figure  11-12.  Gathered  Traces  of  All  Shots  Recorded  at  S2  (top)  and 
SI  (bottom),  Horizontal  Instruments 
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Figure  11-14.  Gathered  Traces  of  E5  (top)  and  E6  (bottom) 
Recorded  at  All  Stations 
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Figure  11-15.  Gathered  Traces  of  E30  (top)  and  E31  (bottom) 
Recorded  at  All  Stations 
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Figure  II- 17.  Gathered  Traces  of  All  Shots  Recorded  at  S31  (top) 
and  S23  (bottom),  Horizontal  Instruments 
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Figure  11-18.  Raw  Traveltimes  of  Phase  I  Secondary  Arrivals 
for  All  Shots  Recorded  at  Each  Station 
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Figure  11-19.  Raw  Traveltimes  of  Phase  I  Secondary  Arrivals 
for  Each  Shot  Recorded  at  All  Stations 
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Figure  11-21.  Raw  Traveltimes  of  Phase  III  Secondary  Arrivals 
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Figure  11-22.  Corrected  Traveltimes  of  Phase  I  Secondary  Arrivals 
for  All  Shots  Recorded  at  Each  Station 


Figure  11-23.  Corrected  Traveltime-  of  Phase  I  Secondary  Arrival 
for  Each  Shot  Recorded  at  All  Stations 
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Figure  11-24.  Corrected  Traveltime s  of  Phase  I  Secondary  Arrivals, 

Horizontal  Instruments 


11-34 


scloneo  ssrvlcss  division 


Figure  11-25.  Corrected  Traveltimes  of  Phase  III  Secondary  Arrivals 
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Figure  11-26.  Reduced  Traveltimes  of  Phase  I  Secondary  Arrivals  for 
All  Shots  Recorded  at  Each  Station 
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Figure  11-27.  Reduced  Traveitimes  of  Phase  I  Secondary  Arrivals 
for  Each  Shot  Recorded  at  All  Stations 
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Figure  11-29.  Reduced  Traveltimes  of  Phase  III  Secondary  Arrivals 


3.  Deconvolution 


In  general,  the  shot  x-ecordings  exhibit  a  reverberatory  charac¬ 
teristic  which  is  fairly  typical  of  data  recorded  in  an  ocean  environment.  This 
effect  is  due  to  the  water  layering  and  tends  to  lengthen  the  refraction  arrivals, 
making  it  difficult  to  detect  separate  refractions.  This  problem  can  be  over¬ 
come  theoretically  by  designing  deconvolution  filters  which  contract  the  arrival's 
duration  by  broadening  its  frequency  band. 

Deconvolution  filters  were  designed  and  applied  to  records  com¬ 
prising  all  shots  recorded  at  one  station  on  the  vertical  and  pressure  channels. 
The  stations  used  were  SI,  S7,  S10,  S23,  and  S3 1.  Figures  11-30  through 
11-34  show  the  input  records,  gate  lengths  used  to  compute  the  autocorrelation 
functions  (and  subsequent  deconvolution  filters),  and  the  deconvolved  records. 

Deconvolution  applied  to  station  S31  is  shown  in  Figure  11-30. 

The  top  record  shows  the  Phase-III  north  shots  aligned  at  a  velocity  of  8.  0 
km/sec.  The  middle  record  shows  deconvolution  results  using  filters  designed 
from  each  trace  over  the  first  gate,  and  the  bottom  record  shows  results  from 
filters  designed  over  the  second  gate.  Broadening  of  the  signal's  frequency 
content  by  the  process  is  apparent.  Water  depth  at  the  station  (1.5  km)  and 
at  the  shots  (3.7  km)  indicates  a  2-way  traveltime  in  the  water  layer  of  about 
2  sec  for  the  station  and  5  sec  for  the  shots.  The  filter  lengths  of  about  0.  5 
sec  were  not  expected  f  "1Jminate  the  bounces  associated  with  the  large  water 
depths  at  most  recording  *  i  ns  and  shot  locations. 
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Figure  11-31  shows  deconvolution  results  using  two  gates  for 
station  S23.  The  water  depth  at  this  station  was  about  0.8  km,  indicating  a 
2-way  traveltime  in  the  water  of  just  over  1  sec.  The  traces  have  been 
whitened  (including  the  noise  preceding  the  signal),  but  little  or  no  improve¬ 
ment  in  phase  identification  has  been  achieved. 

Deconvolution  results  for  station  SI  are  shown  in  Figure  11-32. 
Water  depth  at  this  station  was  relatively  shallow  (0.  34  km),  and  the  noise 
level  was  high.  Whitening  of  the  traces  was  apparent;  however,  little  improve¬ 
ment  in  the  capability  to  identify  masked  secondaries  was  achieved. 

Deconvolution  results  for  stations  S7  and  S10  are  shown  in 
Figures  11-33  and  11-34.  These  show  the  water  bounce  in  deep  water  quite 
clearly.  In  addition  to  using  the  two  short  gates,  shown  in  (b)  and  (c),  a  long 
gate  (17.6  sec)  was  used  to  design  a  longer  deconvolution  filter  (1.73  sec).  The 
result  is  shown  in  (d)  where  the  traces  have  been  shifted  to  vertically  align  the 
initial  P-wave.  Water  depth  at  S7  is  7.  02  km  aud  is  4.  92  km  at  S10,  giving 
a  2-way  traveltime  in  the  water  of  9.  3  and  6. 6  sec,  respectively.  No  water- 
bounce  arrivals  are  possible  in  the  5-  or  6-sec  interval  following  the  initial 
P-wave;  therefore,  clear  secondary  arrivals  would  be  expected  in  that  interval. 
Fragmentary  lineups  are  noted;  however,  they  cannot  be  confidently  picked 
across  the  records  either  before  or  after  deconvolution. 

The  deconvolution  technique  generally  does  not  appear  to 
improve  the  capability  to  pick  secondary  arrivals,  suggesting  that  secondary 
arrivals  were  usually  weak.  This  observation  again  stresses  the  need  to 
design  shot  and  seismometer  arrays  which  allow  the  use  of  first-arrival 
data  to  map  all  refractors. 
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Figure  11-33.  Deconvolution,  Station  S7 
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B.  REFRACTION  ANALYSIS 

The  refraction  analysis  was  fully  described  in  the  Preliminary 
Analysis  Report.  A  brief  description  of  that  analysis  and  results  are  presented 
here  for  comparison  with  other  methods. 

The  appendix  gives  the  statistical  confidence  levels  for  the  param¬ 
eters  of  the  refraction  crustal  model. 

1.  Iterations 

Traveltimes  of  each  first  arrival  were  determined  and  input  to 
a  series  of  computer  programs  which  reduced  the  data  using  standard  refraction- 
analysis  methods.  The  first  program  corrected  the  data  to  a  reference  plane; 
this  correction  was  necessary  because  of  the  large  lateral  velocity  variations 
in  the  uppermost  4  to  5  km  below  sea  level  because  of  large  changes  in  water 
depth  and  structural  environment.  The  reference  plane  was  chosen  to  approx¬ 
imate  the  water  bottom  on  either  side  of  the  Aleutian  Ridge.  Traveltime  plots 
of  the  raw,  corrected,  and  reduced  times  were  generated  by  the  second  program. 
The  third  calculated  least- square  straight-line  fits  to  the  data.  The  final  pro¬ 
gram  resolved  reverse-profile  information  to  obtain  the  velocity,  dip,  and 
depth  of  the  refractors. 

The  data  analysis  was  iterative  in  nature.  Using  information 
from  literature,  the  arrival-time  data  were  reduced  and  crustal  parameters 
estimated.  The  data  were  then  recycled  using  the  output  of  the  previous  iter¬ 
ation  and  were  cycled  through  the  computer  programs  three  times. 

2.  Refraction  Model 

Model  computation  assumes  plane  constant- velocity  dipping 
layers.  Reverse-profile  solutions  were  obtained  from  the  equations  given  in 
Figure  11-35.  Dips  and  depths  of  the  refractors  were  computed  from  measured 
apparent  velocities  and  intercept  times  using  the  output  of  the  final  iteration. 

Line  fits  and  points  used  are  shown  in  Figure  11-36  through  11-39. 
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where 

R  =  horizontal  distance  between  A  and  B 


Vj  =  velocity  in  layer  1 
=  velocity  in  layer  2 

8  -  critical  angle  defined  by  sin  8  -  Vj/V^ 

=  apparent  velocity  in  the  A-to-B  direction  for  waves 
traveling  in  layer  2 

=  apparent  velocity  in  the  B-to-A  direction  for  waves 
traveling  in  layer  2 


t^  =  zero-distance  intercept  time  at  A 

t  =  zero-distance  intercept  time  at  B 

6  =  dip  of  the  interface  in  degrees 

h^  =  vertical  depth  to  the  interface  at  A 

h  =  vertical  depth  to  the  interface  at  B 

ou  r 


Figure  11-35.  Equations  for  Solving  Reverse  Profile 
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Figure  U-37.  Fitted  Single  Profile  of  Phase  I  Shots  Recorded  at  S10 
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The  refraction  model  is  shown  in  Figure  11-40.  Areas  of 
single  (1-way)  coverage  and  double  (reverse)  coverage  are  indicated  on  the 
figure.  The  resolved  Moho  velocity  was  8. 03  km/sec  under  the  Phase  I 
array  and  8.  09  km/sec  under  the  Phase  III  array.  A  reverse-time  mistie 
of  the  Phase  III  reverse  profile  indicated  that  the  segment  under  the  Phase  III 
array  could  not  be  represented  accurately  by  a  single  plane-dipping  layer.  The 
1-way  coverage  under  the  two  shot  arrays  was  resolved  by  assuming  that  the 
Moho  refractor  velocity  was  the  same  as  that  determined  under  the  stations 
for  Phases  I  and  III,  respectively.  Moho  depth  at  the  northern  end  of  the 
Petrel  Bank  was  about  15  km,  which  agrees  with  Shor's  results  for  the  Aleutian 
Basin.  Near  Amchitka,  the  depth  increases  to  over  40  km.  At  the  south¬ 
ern  end  of  the  Phase  I  array,  the  Moho  depth  is  about  12  km,  which  is  con¬ 
sistent  with  normal  oceanic  structure.  In  addition,  two  Pg  refractors  having 
velocities  of  4.  9  and  6.  2  km/ sec  and  depths  of  about  1  and  10  km  were  identified 
on  the  Aleutian  Ridge. 

C.  WAVEFRONT  ANALYSIS 
1.  Wavefront  Methods 

A  wavefront  is  defined  as  the  surface  passing  through  the  most 
advanced  positions  reached  by  a  specific  disturbance  at  any  particular  time. 

A  "wavefront  diagram"  is  therefore  defined  as  a  composite  figure  showing 
the  intersections  of  the  various  wavefronts  with  a  specified  plane  at  succes¬ 
sively  equal  time  intervals.^  The  entire  construction  of  the  wavefront  diagram 
is  based  on  Huygen's  principle.  Applying  this  principle  in  reverse,  the  wave- 
fronts  are  traced  backward  from  the  receiver  to  the  shot.  Using  this  approach 
for  two  in-line  explosions,  one  on  each  end  of  the  line  of  receivers,  wavefronts 
are  determined  for  equal  time  intervals.  The  wavefronts  from  opposing  ex¬ 
plosions,  with  wavefront  times  adding  to  the  reciprocal  time  (RT),  intersect 
on  or  just  below  the  refractor  interface. 
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Plotting  arrival  times  on  a  time-distance  plot  (Figure  11-41) 
is  the  initial  step  in  applying  the  wavefront  method.  By  plotting  the  arrival 
times  from  two  explosions  at  opposite  ends  recorded  by  several  receivers 
in  a  line,  a  reverse  time  can  be  determined.  Equal  traveltimes  over  the 
whole  spread,  from  one  shotpoint  to  the  other,  are  necessary  for  the  method 
to  produce  a  reliable  model;  i.e.  ,  the  reverse  traveltimes  must  agree.  Reli¬ 
ability  of  the  reverse-time  data  determines  the  reliability  of  he  path  inter¬ 
pretation  and  the  refractor  velocity  determinations. 

After  the  time-distance  curves  are  plotted,  x-positions 
corresponding  to  regular  increments  of  arrival  time  at  some  convenient  At 
are  marked  off  on  the  curves.  The  x-positions  are  then  transferred  to  the 
surface  line  of  the  receivers.  The  reversed  curve  is  marked  off  in  similar 
increments  of  time;  however,  these  are  at  values  of  reverse  time  minus  the 
integral  time  values.  From  this  point,  different  variations  of  the  wavefront 

method  may  be  applied.  The  two  methods  used  were  the  general  wavefront 

ft  9 

method  and  wavefront  targeting. 

Using  the  general  wavefront  method,  arcs  of  small  circles  of 
radii  equal  to  At  •  V 2At  ■  Vj,  3At  •  Vj,  etc.  ,  were  drawn  from  the  surface 
position  of  the  arriving  wavefronts  back  toward  the  event.  Tangents  to  these 
sets  of  wavelets  represent  the  wavefront  system  emanating  from  one  of  the 
explosions.  This  same  procedure  was  used  to  produce  a  wavefront  system 
from  the  other  explosion.  Since  the  time  values  of  one  set  of  wavefronts  are 
chosen  at  the  regular  increments  of  At  and  the  other  set  of  wavefronts  have 
time  values  of  reverse  time  minus  n  •  At,  the  wavefronts  intersect  at  some 
depth  so  that  the  sum  of  the  times  of  intersecting  wavefronts  represents  a  total 
time  equal  to  the  reverse  time.  This  point  of  intersection  is  a  point  on  the 
refractor.  The  horizontal  diagonal  formed  by  the  intersection  of  consecutive 
wavefronts  is  proportional  to  the  true  velocity  of  the  refractor.  This  velocity 
equals  the  measured  length  of  this  diagonal  divided  by  the  At  used. 

The  second  wavefront  method  used  was  wavefront  targeting. 

This  method  permits  the  rapid  construction  of  the  refractor  from  reversed 


II  56 


solano*  s*rvloos  division 


data.  The  method  consists  of  reproducing  two  complete  complementary  wave- 
front  systems  for  a  given  refractor,  including  image  portions  below  the 
refractor  and  imaginary  portions  above  the  surface.  These  systems,  as  in  the 
general  wavefront  method,  intersect  on  or  slightly  below  the  refractor.  One 
wavefront  is  constructed  by  drawing  arcs  of  circles  with  centers  at  positions 
of  arrival  time  and  with  radii  equal  to  the  difference  between  the  arrival  time 
and  the  wavefront  time  multiplied  by  Vj.  The  wavefront  is  formed  by  connec¬ 
ting  the  arcs  by  tangents  to  the  arcs.  The  other  wavefront,  with  time  equal  to 
RT  -  time  of  the  other  wavefront,  is  constructed  in  a  similar  manner  using 
arrival  times  from  the  other  source  as  the  centers  for  the  arcs. 

A  wavefront  target,  concentric  circles  of  radii  equal  to  At  •  Vp 
2At  •  V},  3At  •  V  i ,  etc.,  was  drawn  on  a  transparent  plastic  sheet.  A  small 
hole  was  punched  in  the  center  of  the  target  to  accommodate  a  sharp  pencil  point. 
If  a  circle  is  tangent  to  both  wavefronts  of  the  complementary  system,  the 
center  of  the  circle  should  be  another  point  on  the  refractor.  The  target 
was  maneuvered  so  that  each  successive  circle  was  tangent  to  both  wavefronts, 
and  the  location  of  the  center  was  marked  each  time.  The  series  of  points 
produced  in  this  manner  were  spaced  at  intervals  of  At  along  the  refracting 
interface.  The  measured  length  between  each  pair  of  successive  points 
divided  by  At  readily  gives  the  refractor  velocity  for  that  interval. 

2.  Wavefront  Model 

Since  the  Aleutian  Islands  Experiment  was  shot  with  the  re¬ 
ceivers  on  the  water  bottom  and  the  explosions  effectively  on  the  surface, 
the  arrival  times  were  corrected  to  a  reference  plane  in  order  to  have  source 
and  receiver  on  the  same  plane.  The  corrected  arrival  times  were  plotted 
as  time  (sec)  vs  distance  (km),  giving  a  time-distance  plot  (Figure  11-41). 

Since  there  was  no  explosion  located  at  either  receiver  location,  the  exact 
reverse  time  could  not  be  measured.  Therefore,  to  determine  the  reverse 
time,  the  line  through  the  plotted  points  on  the  time-distance  curve  must  be 
projected  to  the  distance  which  locates  the  source;  the  corresponding  time 
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represents  the  reverse  time.  Phase  I  data  were  used  (Figure  11-41).  The 
reverse  time  for  E34  to  E7  equals  33  sec. 


Data  recorded  at  six  locations  (S5,  S6,  S7 ,  S8,  S9,  and  S10) 
were  used.  First  arrivals  from  E7  to  SI  and  S2  were  refracted  by  a  shallow 
interface;  therefore,  SI  and  S2  cannot  be  used,  since  only  one  refractor  at  a 
time  can  be  mapped  using  these  methods.  The  arrival  time  from  E34  to  S10 
cannot  be  used  because  it  was  a  reflection  (or  refraction  from  a  shallow 
interface)  and  only  refractions  from  the  same  interface  can  be  used  in  this 
method.  Although  only  a  few  arrival  times  from  the  same  interface  were 
available,  the  plotted  points  were  connected,  allowing  points  to  be  picked  at 
a  convenient  At,  say  1  sec;  those  points  were  spotted  at  their  corresponding 
distances  along  the  surface  line.  Thus,  the  number  of  times  available  were 
increased,  giving  better  coverage. 

After  the  times  were  transferred  to  the  surface  line  (in  this 
case,  the  reference  plane),  arcs  were  drawn  with  radii  equal  to  1  sec  •  Vp 
2  sec  .  Vp  3  sec  •  V^,  etc.  ,  where  =  6.  4  km/sec.  In  places,  arcs  drawn 
from  adjacent  points  overlapped,  thus  making  drawing  tangents  to  these  arcs 
impossible.  This  problem  occurs  in  the  area  of  S8  and  S9.  In  order  to  obtain 
a  wavefront  solution,  the  curve  connecting  the  arrival  times  on  the  time-distance 
plot  was  smoothed,  thereby  adjusting  the  estimated  arrival  times  and  allowing 
arcs  to  be  properly  drawn. 

By  employing  the  general  wavefront  method,  wavefronts  from 
E7  with  times  of  16  sec,  17  sec,  18  sec,  ....  25  sec,  26  sec  were  used. 

These  wavefronts  intersected  with  wavefronts  from  E34  with  times  of  17  sec, 

16  sec,  15  sec,  ....  8  sec,  7  sec  to  form  points  on  the  refractor.  This  pro¬ 
cedure  produced  the  model  shown  in  Figure  11-42  with  a  refractor  velocity 
(V^)  equal  to  7. 88  km/sec. 

This  model  slopes  downward  toward  E7  with  an  average  dip 
of  13°.  The  shallowest  point  was  under  S10  (13.5  km),  and  the  deepest  point 
was  under  S5  (34  km);  all  depths  were  measured  from  the  4.  5-km  reference 
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plane.  The  portion  of  the  refractor  determined  by  the  intersection  of  wave- 
fronts  with  times  to  the  north  of  S5  is  projected  points  and  therefore  was 
not  necessarily  reliable.  This  includes  the  last  two  segments  on  the  north  end 
of  the  refractor.  The  last  three  segments  on  the  south  end  also  were  not 
necessarily  reliable,  since  the  wavefronts  for  arrival  times  of  27  sec,  28  sec, 
and  29  sec  were  based  on  only  1-way  data.  (Arrival  times  from  E34  were 
not  available  in  this  area,  so  there  were  no  wavefronts  to  intersect  with  those 
from  E7).  In  order  to  get  the  points  south  of  S9,  a  radius  equal  to  At  times 
the  already  determined  refractor  velocity  drawn  with  the  center  at  the  extreme 
point  on  the  refractor  intersects  with  1-way  wavefronts  from  E7.  In  this  man¬ 
ner,  the  refractor  can  be  extended  to  S10.  The  model,  therefore,  is  most 
reliable  in  the  area  between  S8  and  a  point  10  km  north  of  S6. 

Using  the  wavefront  targeting  method,  a  wavefront  system, 
which  included  a  wavefront  from  E34  with  a  time  of  12  sec  and  a  wavefront 
with  a  time  of  21  sec  from  E7,  was  developed.  By  maneuvering  the  wavefront 
target  so  that  each  successive  circle  was  tangent  to  both  complements  of  the 
system,  the  wavefront  targeting  model  was  produced  (Figure  11-43).  This 
method  gave  a  refractor  velocity  of  7.  98  km/sec. 

This  model  was  very  similar  to  that  produced  by  the  other 
method.  The  dip  trend  was  in  the  same  direction  as  that  of  the  previous  model; 
the  average  dip  was  12°,  which  agrees  with  refractor  dip  determined  by  per¬ 
turbation  techniques  over  the  same  segment.  The  shallowest  point  (13.5  km) 
was  3.  5  km  south  of  S9,  and  the  deepest  point  (29.  5  km)  was  7  km  north  of  S6. 

Both  of  the  wavefront  methods  produce  models  of  the  refractor 
which  are  reliable;  however,  to  achieve  the  best  model  possible,  a  more  ac¬ 
curate  method  of  measuring  the  reverse  traveltime  must  be  used.  Ideally, 
there  should  be  a  receiver  at  each  source  location  to  directly  record  the 
entire  traveltime  from  both  directions.  Also,  there  should  be  more  receivers 
along  the  line  to  directly  record  more  arrival  times.  These  would  give  better 
coverage  over  the  spread  and  thus  give  better  reliability  to  the  interpretation. 
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Figure  11-41.  Reverse  Profile  Plot  Between  E34  and  E7 ,  Phase 
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Figure  11-42.  Wavefront  Model  Using  General  Wavefront  Method 


D.  TIME-TERM  ANALYSIS 


1.  Time-Term  Theory 

The  basic  equation  of  the  time-term  method  of  seismic  refraction 
analysis  is  derived  from  the  geometry  of  a  refraction  model  consisting  of  two 
layers  separated  by  a  plane  interface.  The  solution  of  a  system  of  such  equa¬ 
tions  for  experimental  traveltime  data  from  a  2-layered  structure  yields  a 
seismic  velocity  in  the  lower  layer  and  the  contour  of  the  interface.  The  con¬ 
tour  is  given  in  terms  of  the  depth  to  the  interface  beneath  each  shot  and  station 
of  the  experimental  array. 

In  Figure  11-44,  A  is  the  distance  between  shot  A  and  station  B; 

H&  and  are  the  perpendiculars  from  the  refracting  layer  to  A  and  B;  0  is 
the  angle  of  incidence  of  the  critically  refracted  ray;  $  is  the  dip  angle  of  the 
refractor;  and  are  the  seismic  velocities  in  the  upper  layer  and  refrac¬ 
tor,  respectively;  and  dp  d2>  and  d-j  are  the  three  segments  of  the  ray  path 
from  A  to  B.  The  traveltime  t  of  the  ray  from  A  to  B  is 


Also  from  Figure  11-44, 

H 
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Figure  11-44.  Ray  Path  Diagram  for  Critically  Refracted  Ray 


Applying  the  relation  Vj/V  =  sin  0  to  Equation  2-5  yields 


t  =  (H  +  H.)  t 

a  b  v  v 

1  v  2 


(2-6) 


It  may  be  shown  that  for  dip  angles  normally  encountered,  the 
error  in  (H&  +  H^)  which  arises  from  setting  cos  $  =  1  in  Equation  2-6  is 
negligibly  small.  Making  this  approximation  in  Equation  2-6  results  in  the 
time-term  equation 


t  =  P  +  q  + 


(2-7) 


where 
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Therefore,  in  the  tirrie-term  equation,  the  time-term  p  for  a  given  shot  and 
the  time-term  q  for  a  given  station  are  characteristic  of  the  velocities  Vj  and 
V 2  and  the  depth  H  of  the  refractor  beneath  the  shot  or  station  site. 

The  refractor  velocity  in  Equation  2-7,  previously  considered 

constant,  may  be  replaced  with  a  refractor  velocity  which  is  a  linear  function 

of  the  distance  A.  Changes  in  V.  (A)  may  be  attributed  to  a  velocity  gradient 

1 1  L 

in  the  refracting  medium.  This  velocity  function  may  be  written  as 


V  (A)  s  V  +  V,  A 
2  a  b 


(2-9) 
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In  an  array  consisting  of  m  shots  and  n  stations,  a  maximum  of 
m  x  n  equations  of  the  form  of  Equation  2-13  may  be  written.  However,  if  no 
shot  and  station  locations  coincide,  this  system  will  not  have  a  unique  solution. 
As  inspection  of  Equation  2-13  shows,  an  arbitrary  constant  a  may  be  subtrac¬ 
ted  from  each  shot  time-term  and  added  to  each  station  time-term  without 
affecting  the  equation  otherwise.  To  overcome  this  difficulty,  one  of  the  time- 
terms  may  be  set  equal  to  0  arbitrarily  and  the  constant  a  reintroduced  after 
the  system  of  equations  has  been  solved. 


The  system  of  Equations  2-13  for  an  array  of  m  shots  and  n 
stations  may  be  represented  in  matrix  form  as 

AB  =  T  -  kL  D  -  k2  E  (2-15) 

The  matrix  B  is  a  column  matrix  of  m  +  n  -1  elements,  the  first  m  elements 

bj  (j  =  1,  2,  ....  m)  being  the  shot  time-terms  and  the  last  n  -  1  elements 

bj  (j  =  m  +  1,  m  +  2,  ....  m  +  n-1)  being  the  station  time-terms.  The  last 

station  time-term  t>m+n  has  been  set  equal  to  0,  as  previously  explained.  The 

matrix  A  is  a  matrix  of  L  rows  and  m+n-1  columns,  where  L  *  (m  x  n)  is 

the  number  of  observations  obtained  from  the  array.  Its  elements  a.,  are 

1J 

either  l's  or  0's,  depending  on  whether  or  not  the  i*h  observation  was  either 
initiated  or  recorded  at  the  j**1  site.  The  matrix  T  is  a  column  matrix  of  the  L 
traveltimes  observed;  D  is  the  corresponding  matrix  of  the  A  for  each  obser¬ 
vation;  and  E  is  the  corresponding  matrix  of  the  squared  A  for  each  observation. 
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where  Q  =  A'A,  P  =  A'T,  R  =  A'D,  and  S  =  A'E.  Multiplying  on  the  left  by 
Q_1,  the  inverse  of  Q,  yields 

B  =  X  -  kj  Y  -  k2  Z  (2-17) 

where  X  =  Q-1P,  Y  =  Q_1R,  and  Z  =  Q-1S. 


Equation  2-17  gives  the  solution  for  the  m  +  n  -  1  time-terms, 
but  the  values  of  and  are  not  yet  determined.  To  find  the  best-fitting 
values  of  k^  and  k^  ,  Equation  2-17  is  first  substituted  into  Equation  2-15, 
yielding 

C  -  kjF  -  k2  G  =  0  (2-18) 


where 


C  =  AX  -  T 
F  =  AY  -  D 
G  =  AZ  -  E 


(2-19) 


The  matrices  C,  F,  and  G  are  each  column  matrices  with  L 
elements.  To  find  the  least-squares  value  of  k^  and  k^,  it  is  necessary  to 
consider  the  equivalent  system  of  L  equations  which  can  be  formed  from 
Equation  2-18.  Considering  the  i^1  such  equation  and  using  6^  to  designate 
the  residual  which  results  from  fitting  experimental  data, 

c.  -  k  f.  -  k  g.  =  6  .  (2-20) 

l  1  l  2  ®i  i 

where  c^,  f^,  and  g^  are  the  i**1  elements  of  the  matrices  C,  F,  and  G,  respec¬ 
tively.  For  k^  and  k^  to  best  fit  the  data  in  the  least-squares  sense,  the  sum 
of  the  squared  residuals  must  be  a  minimum  with  respect  to  k^  ? nd  k^. 

The  expression  for  the  sum  I  of  the  squared  residuals  is 

L  L  2 

1  *  .£  sf  ■  .£  [ci  -  <Vi  +  k2*i)]  (2-21) 
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Taking  the  partial  derivative  of  I  with  respect  to  and  equating  the  result  to 
0  yields  the  equation 


k.E  fi2  +  k2  z  E 

i  =  1  i  =  1  i  =  1 


c.  f. 
1  1 


(2-22) 


or,  employing  the  summation  convention  for  repeated  subscripts, 


k.  f.  f.  +  k  f.  g.  =  c.  f. 
lii  2  i  ®i  ii 


(2-23) 


Similarly,  taking  the  partial  derivative  of  I  with  respect  to  k^  and  equating 
the  result  to  0  yields 


k  f.  g.  +  k  g.  g.  =  c.  g. 
lii  2  i  l  i  i 


(2-24) 


Equations  2-23  and  2-24  may  be  solved  for  k^  and  k^,  yielding 


kl  = 


(ci  ^i)  (Si  8i)  ”  (ci  Si )  (fj  8i) 

<Vi>  <gigi>  '  <£i  *i»2 


(2-25) 


and 


V 


(f.  f.)  (c.  g.)  -  (f.  g.)  (c.  f.) 
i  i  '  i  *i  '  i  ®i  '  i  i 


<Vi»  <gigi’  -  <fi  gt>‘ 


(2-26) 


Equations  2-25  and  2-26  give  the  least-squares  values  for  k^ 
and  k^  in  terms  of  known  quantities.  These  values  may  then  be  substituted  into 
Equation  2-17  to  solve  for  the  m  +  n  -  1  time-terms  and  into  Equation  2-14  to 
determine  the  two  parameters  of  the  refractor  velocity. 
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2.  Time-Term  Model 


The  Moho  refraction  arrivals  from  Phase  I,  Phase  III,  and 
four  land  stations  were  processed  through  the  time-term  computer  program. 
Input  to  the  program  were  the  time-distance  pairs  for  each  shot- station  ob¬ 
servation  shown  in  Table  II - 3 .  The  time-term  for  S8  was  arbitrarily  set  equal 
to  0,  thereby  reducing  the  order  of  the  matrix  by  1  and  removing  the  singularity 
in  the  coefficient  matrix.  After  solving  the  system  of  equations,  the  normal¬ 
izing  constant  a  was  reintroduced  by  assuming  that  shot  E24  and  station  S2 
were  sufficiently  close  together  to  have  identical  time-terms  (the  separation 
was  2.9  km).  The  time-term  results  can  be  presented  in  either  of  two  ways: 
shown  in  sec  or  converted  to  depths  (H),  where 


In  this  report,  time-terms  are  converted  to  depths  using  a 
crustal  velocity  (V^)  of  6.4  km/sec  for  comparison  with  models  developed 
using  other  techniques  (i.e.  ,  refraction,  wavefront,  and  model  perturbation). 

The  time-term-model  solution,  Moho  velocity,  and  rms  resid¬ 
ual  (for  all  observations)  are  shown  in  Figure  11-45.  An  examination  of 
these  results  shows  (in  addition  to  observational  scatter)  that  the  refractor 
is  dipping  toward  Amchitka  and  that  the  deepest  part  of  the  structure  lies 
under  the  Phase  III  stations.  This  is  in  general  agreement  with  the  refraction- 
nodel  and  perturbation-model  results;  however,  a  large  vertical  mistie 
beneath  Amchitka  (S40)  and  an  apparent  lack  of  continuity  between  Phase  III 
north  shots  and  stations  were  noted. 
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SHOT/STATION  EXPERIMENTAL  DATA 
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(X4  vertical  exaggeration) 


The  OBS  unit  and  shot  arrays  were  deployed  in-line,  while  the 
land  stations  atAdak,  Atka,  Shemya,  and  Amchitka  were  generally  off-line 
(Amchitka  being  most  nearly  in-line  with  the  profile).  Since  scatter  in  some 
of  the  land-station  observations  was  noted,  the  OBS  data  were  processed  in¬ 
dependently;  however,  the  Phase  I  and  Phase  III  operations  were  conducted  at 
different  times  which  produced  a  singular  observational  matrix.  The  USC&GS 
station  at  Adak  had  very  consistent  arrival  times  and  was  included  to  provide 
a  connection  between  Phase  I  and  Phase  III  observations.  The  resulting 
time-term-model  solution,  Moho  velocity,  and  rms  residual  (for  all  obser¬ 
vations)  are  shown  in  Figure  11-46.  The  results  were  not  significantly  different 
from  previous  results  (Figure  11-45). 

Phase  I  and  Phase  III  data  were  processed  independently  through 
the  time-term  program;  and  although  better  continuity  could  be  obtained  between 
Phase  III  north  shots  and  stations  by  adjusting  the  Phase  III  constant  (a),  the 
structure  shown  under  Amchitka  was  not  continuous  and  exhibited  essentially 
the  same  structural  pattern  as  before 

The  constrained  velocity  solutions  for  each  of  these  cases 
(7.  85  through  8.  20  km/sec  in  increments  of  0.  5  km/sec)  produced  some 
tilting  of  the  structure  under  the  shots  with  respect  to  structure  under  the 
stations.  However,  the  same  basic  structural  pattern  was  observed  (i.  e.  , 
vertical  mistie  beneath  Amchitka  and  poor  continuity  of  structure  beneath 
Phase  III  north  shots  and  stations). 

3.  Theoretical  Example 

To  examine  the  structure  yielded  by  the  time-term  method  for 

a  simplified  version  of  the  final  crustal  model,  the  2 -layer  crustal  model  of 
Figure  11-47  (the  interface  of  which  is  indicated  by  the  solid  line)  was  input  to 
the  theoretical  traveltime  program.  The  resulting  theoretical  traveltimes 
were  then  input  to  the  time-term  program,  yielding  the  time-term  depths 
indicated  by  the  points  in  Figure  11-47. 
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Figure  11-46.  Time-Term  Model  Using  Data  from  Phase  I,  Phase  III, 
and  Adak  Land  Station  (X4  vertical  exaggeration) 
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Figure  11-47.  Time-Term  Results  Using  Theoretically  Computed  Traveltimes. 

The  interface  of  the  theoretical  model  (with  dip  angles  less  than 
11.8°)  is  indicated  by  the  solid  line  (X4  vertical  exaggeration). 


Table  II-4  shows  the  shot/station  pairs  used  in  computing  the 
time-term  depths.  Inspection  of  Figure  11-47  and  Table  II-4  shows  that 
generally  poor  results  are  obtained  by  the  time -term  method  for  those  stations 
under  the  root  of  the  model,  i.  e  ,  stations  S23,  S25,  S28,  S30,  and  S31  and 
those  shots  which  connect  with  them.  Relatively  good  agreement  between 
the  true  model  and  the  time- term  depths  is  obtained  for  those  sites  not 
associated  with  the  deepest  parts  of  the  model  structure,  i.  e.  ,  those  sites 
to  the  left  of  E26. 

Since  some  of  the  disagreement  between  the  true  and  the  time- 

term  depths  might  be  attributable  to  violation  of  the  assumption  that  the  cosines 

of  the  dip  angles  may  be  approximated  by  unity,  another  model,  with  shallower 

dip  angles,  was  input  to  the  time-term  program.  The  model  and  the  time- 

term  depths  derived  from  it  are  shown  in  Figure  11-48.  The  largest  dip  angle 

for  the  model  of  Figure  11-47  was  11.8°;  the  largest  dip  angle  for  the  model 

of  Figure  11-48  Wo.s  5.2°.  Although  the  time-term  fit  is  somewhat  better  in 

this  example,  substantial  deviation  from  the  true  model  still  occurs  for  the 

stations  located  over  the  root  of  the  structure  and  the  shots  associated  with 

12 

them.  This  deviation  supports  the  conclusion  of  previous  investigations 
that  the  time-term  method  is  not  fully  suitable  for  deriving  the  structure  under 
a  site  when  the  interface  within  the  cone  of  critically  refracted  rays  from  that 
site  deviates  substantially  from  planarity. 

E.  TRAVELTIME  ANALYSIS 

1,  Traveltime  Program 

A  traveltime  program  was  developed  to  compute  theoretical 
traveltimes  for  crustal  models  with  layers  which  change  dip.  The  program 
contains  a  plot  option  so  that  the  model  may  be  plotted  by  computer,  if  de¬ 
sired;  a  computer -generated  plot  of  a  typical  model  is  shown  in  Figure  11-49. 
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Figure  11-48.  Time-Term  Results  Using  Theoretically  Computed  Traveltimes. 

The  interface  of  the  theoretical  model  (with  dip  angles  less  than 
5.2°)  is  indicated  by  the  solid  line  (X4  vertical  exaggeration) 
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Figure  11-49.  Example  of  Output  from  Program  Which  Computes  Traveltimes 
for  Crustal  Models  with  Layers  that  Change  Dip 


The  model  is  input  to  the  program  by  describing  the  interface  above  each  layer 
in  terms  of  the  horizontal  distance  and  vertical  depth  from  the  origin  (0  in 
Figure  11-49)  to  each  point  on  the  interface  at  which  a  change  of  dip  angle 
occurs  (e.  g.  ,  the  point  B  on  the  interface  above  layer  IV,  where  the  dis¬ 
tance  and  depth  are  OP  and  PB,  respectively).  The  program  will  accept 
a  total  of  nine  layers  and  a  half-space  and  a  maximum  of  80  breakpoints 
(such  as  point  B)  in  any  one  interface.  Also,  a  P-wave  velocity  for  each 
layer  is  read  in. 

The  locations  of  both  sources  (e.  g.  ,  points  S)  and  receivers 
(e.  g.  ,  points  R)  are  specified  in  terms  of  their  horizontal  distances  from  the 
origin  and  the  number  of  the  interface  on  which  they  lie;  sources  and  receivers 
may  be  specified  to  lie  on  separate  interfaces.  One  hundred  source  locations 
and  100  receiver  locations  may  be  specified. 

Ray  paths  are  computed  by  geometrical  acoustics;  rays  are 
assumed  to  enter  and  exit  the  refracting  layer  at  the  critical  angle,  following 
the  contours  of  the  interface  above  the  refractor  between  entry  and  exit.  Ray 
paths  from  the  source  S  to  the  refractor  (layer  V)  and  from  the  refractor  to 
the  receivers  R  are  also  shown  in  Figure  11-49.  For  a  given  input  model,  in 
addition  to  the  options  described  above,  the  program  may  be  recycled  to 
compute  traveltimes  for  additional  source  and  receiver  locations  or  for 
different  refracting  layers. 

2.  Perturbation  Model 

The  traveltime  program  was  used  to  perturb  the  crustal  model 
across  Amchitka  to  give  a  best  fit  to  the  observed  traveltimes.  The  technique 
follows. 
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Given  an  estimate  of  the  2 -dimensional  model,  the  theoretical 
traveltimes  for  each  shot/ station  pair  are  computed.  In  addition,  the  program 
computes  a  residual  (experimental  traveltime  minus  theoretical  traveltime)  for 
each  shot/station  pair  having  an  experimental  observation  and  subsequently 
calculates  the  mean-absolute- value  and  the  root-mean- square  of  the  residuals 
to  give  a  measure  of  how  well  the  theoretical  model  fits  the  experimental 
traveltimes.  An  inspection  of  the  traveltime  residuals  allows  determination 
of  areas  along  the  refracting  interface  where  the  structure  should  be  flexed  to 
give  a  closer  fit  to  the  observed  traveltimes.  The  method,  of  course,  considers 
the  actual  ray  paths  involved.  The  flexed  model  is  input  to  the  computer; 
theoretical  traveltimes  are  determined;  and  residuals  are  computed  as  before. 
The  model  is  then  flexed  again  and  the  process  repeated  as  necessary  to  pro¬ 
duce  a  model  which  fits  the  experimental  observations. 

The  initial  model  was  the  refraction  model,  with  the  refracting 
interface  connected  across  areas  where  no  coverage  was  obtained  in  the  stan¬ 
dard  refraction  analysis.  Average  crustal  velocity  was  taken  as  6.4  km/sec. 
Shown  in  Figures  11-50  and  11-51  are  the  initial  model  and  the  series  of  model 
perturbations.  Mean-absolute-values  (mav)  and  root-mean-squares  (rms)  of 
the  traveltime  residuals  for  each  model  are  given  in  Table  II -  5 .  The  model 
was  perturbed  eight  times. 

Due  to  the  contour  of  the  refracting  interface  (which  is  com¬ 
prised  of  many  straight-line  segments),  a  theoretical  observation  was  not 
obtained  at  shot  E8  for  iterations  3,  4,  and  5,  which  excluded  the  exper¬ 
imental  observations  of  E8  to  S2,  S5,  S6,  S7,  S8,  S9,  and  S10  from  those  mav 
and  rms  calculations.  This  resulted  in  an  increase  of  those  values  on  itera¬ 
tion  6  for  which  a  theoretical  observation  was  obtained  at  shot  E8. 
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Figure  11-51  (iteration  8)  shows  the  final  model,  where  the 
mav  was  0.45  sec  and  the  rms  was  0.  56  sec.  Additional  model  perturbations 
could  undoubtedly  improve  these  values  somewhat;  however,  it  is  felt  that 
they  essentially  represent  the  contributions  due  to  observational  scatter  in 
the  experimental  data.  The  model-perturbation  method  of  crustal  determination 
gives  a  good  representation  of  the  contour  of  the  refractor  across  the  profile, 
is  not  restricted  by  various  assumptions  necessary  to  implement  other  methods 
(e.  g.  ,  standard  refraction  and  time-term),  and  allows  a  direct  measure  of  how 
well  the  model  fits  the  observations.  The  final  perturbation  model  (iteration  8) 
is  a  2-dimensional  determination  from  in-line  profile  recordings.  Although 
the  crustal  parameters  cannot  be  uniquely  determined  from  the  available  data, 
the  model  gives  a  good  fit  to  the  seismic  observations.  The  deepest  part  of  the 
structure  lies  under  the  Petrel  Bank,  and  the  root  of  the  structure  is  located 
60  km  north  of  Amchitka  under  Semisopochnoi  Island. 


Table  II- 5 

TRAVELTIME  RESIDUAL  VALUES 


mav 

(sec) 

rms 

(sec) 

Refraction  Model 

0.61 

0.77 

Iteration  1 

0.63 

0.78 

Iteration  2 

0.  53 

0.65 

Iteration  3 

0.  50 

0.62 

Iteration  4 

0.  43 

0.54 

Iteration  5 

0.  42 

0.53 

Iteration  6 

0.  48 

0.60 

Iteration  7 

0.46 

0.59 

Iteration  8 

0.45 

0.  56 
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Figure  11-50.  Model-Perturbation  Results  Showing  Initial  (refraction) 
Model  and  Iterations  1  Through  4 
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SECTION  III 


EFFECT  OF  CRUSTAL  STRUCTURE  ON  REGIONAL  TRAVELTIMES 

This  section  presents  results  of  calculations  made  to  determine 
the  effect  of  the  Amchitka  crustal  structure  on  traveltimes  observed  at  regional 
distances. 

Utilizing  the  computer  program  described  in  Section  II,  sub¬ 
section  E,  traveltime  curves  were  determined  for  the  symmetrical  model 
shown  in  Figure  III- 1.  The  model  consists  of  a  standard  (JB)  33-km  crust 
with  an  average  crustal  velocity  of  6.05  km/sec  overlying  a  simplified  upper 
mantle  structure  made  up  of  several  concentric  constant- velocity  layers.  Ray 
paths  are  shown  for  the  refraction  in  layer  7,  which  is  the  first  arrival  in  the 
20°  distance  range.  The  calculated  traveltime  curve  for  this  refraction  is 
shown  in  Figure  III  -  2 . 

Next,  the  model  was  modified  simply  by  replacing  the  standard 
crust  on  the  left  side  of  the  model  with  the  Amchitka  Crust  (the  upper  mantle 
structure  being  unchanged)  to  give  the  dual  model  shown  in  Figure  III  - 1 . 
Traveltime  curves  were  computed  for  the  same  refracting  layer;  these  curves 
are  presented  in  Figure  III  -  2  for  comparison  with  the  symmetrical  model 
results.  The  shape  of  the  traveltime  curve  essentially  reflects  the  inverted 
contour  of  the  Moho  refractor  across  Amchitka.  In  addition,  the  traveltime 
curve  is  0.  5  to  2  sec  earlier  than  that  for  the  model  with  the  standard  crust, 
due  primarily  to  differences  in  the  crustal  velocity  of  the  Amchitka  model 
(6.4  km/sec)  and  the  JB  model  (6.05  km/sec). 

For  an  event  20°  from  Amchitka  recorded  along  the  OBS  pro¬ 
file  (or  for  a  series  of  events  along  the  OBS  profile  recorded  at  20°  distance), 
a  traveltime  variation  of  about  1.  5  sec  is  observed  because  of  changes  in  the 
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crustal  structure  along  that  profile.  Given  the  true  upper  mantle  structure 
of  the  region  (and  resultant  traveltime  curve),  variations  of  ±  0.7  sec  are 
observed  if  changes  in  crustal  structure  are  neglected.  These  variations 
are  quite  often  interpreted  as  observational  scatter.  While  these  variations 
may  be  small  compared  to  traveltime  anomalies  resulting  from  incorrect 
upper  mantle  models,  they  are  not  insignificant. 


I 


t 


Figure  III- 1 .  Crustal  Refraction  Ray  Paths  for  Symmetrical  Model,  Standard 
Crust  (top),  and  Dual  Model,  Amchitka  Crust  with  Standard 
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Figure  III  -  2 .  Reduced  Traveltimes  for  Symmetrical  Model  and  Dual  Model 
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A  study  was  undertaken  to  arrive  at  statistical  confidence 
limits  for  the  parameters  of  the  determined  refraction  crustal  model. 

Given  the  interval  (a,  b)  denoted  by  the  limits  a  and  b  and  a  confidence  level 
N,  confidence  limits  stated  for  a  quantity  may  be  interpreted  as  follows:  the 
interval  (a,  b)  will  cover  the  desired  quantity  100  N-percent  of  the  time. 

For  example,  suppose  it  is  desired  to  place  0.95-percent  con¬ 
fidence  limits  on  the  mean  of  a  normal  distribution,  given  four  independent 
samples  from  that  distribution  having  a  variance  equal  to  9.  Using  the  sam¬ 
ple  mean  of  the  observation  X  as  the  estimate  of  the  true  mean  |J  (this  is  the 
maximum -likelihood  estimate  of  ^i),  it  can  be  shown  that 

p|  X  -  2.94  s  u  SX  +  2.94}=  0.95  (A-  1) 

If  the  four  samples  taken  are  such  that  X  =  2.7,  then  a  =  -0.  24  and  b  =  5. 64 
so  that  (-0.24,  5.64)  is  the  95-percent  confidence  interval  on(j.  For  a 
different  set  of  four  samples,  the  values  of  a  and  b  are  different. 

The  first  step  in  arriving  at  confidence  limits  for  the  param¬ 
eters  of  the  crustal  model  was  the  fitting  of  the  points  of  the  time-distance 
plots  to  straight-line  segments.  Ordinarily,  this  is  done  by  a  simple  least- 

squares  fit,  assuming  that  the  distance  errors  in  the  data  are  much  smaller 

13 

than  the  time  error.  Since  this  condition  cannot  be  assumed  to  hold  for  the 
present  data,  it  was  decided  to  employ  a  line  fit  which  assumes  both  distance 
and  time  data  to  be  subject  to  error  and  to  compare  the  results  with  standard 
least-squares-error  fits. 

The  equation  for  a  straight  line  has  the  form 

t  =  7+  b  (X  -  X)  (A-2) 


A-  1 


soli 


irvil 


livisli 


It  is  seen  from  Equations  A-6,  A-7,  and  A-8  that  it  is  not 
2  2 

necessary  to  know  both  a.  and  a  (the  standard  deviations  of  <X.\  and  it.l  ) 

4  x  2  \  l  lf 

but  only  their  ratio.  Given  the  result  of  Equation  A-8,  can  be  estimated 
from  the  N  data  points  by 


N 

S  (xi  -  x)2 

i  =  l 


N 


i=  1 


(X.  -  X)  (t. 


t) 


N  -  2 


(N  -  2)b 


(A-9) 


Using  Equation  2-2,  it  is  seen  that  tQ  is  given  by 

t  =  t  -  b  X  (A- 10) 

when  the  line  is  required  to  pass  through  the  mean  of  the  observations. 

Therefore,  the  line  which  best  fits  the  data  when  both  variables  are  in  error 
has  a  slope  given  by  Equation  A-8  and  an  intercept  given  by  Equation  A-  10. 

For  the  problem  under  consideration,  X  corresponds  to  the  dis¬ 
tance  between  signal  source  and  receiver  and  t  to  the  measurement  of  epoch. 

It  should  be  noted  that  the  distance  error  arises  from  three  considerations: 
curvature  of  the  earth,  error  in  longitude,  and  error  in  latitude. 

The  circumference  of  the  earth  is  21,600  nautical  miles.  If  a 
great  circle  arc  is  used,  the  plane  angle  subtended  for  a  400-nautical  mile  arc 
is  0.  1 16  radians  or  6.  65°  .  The  chord  length  is  then  within  0.  5  percent  of  the 
arc  length.  Hence,  it  will  be  assumed  that  the  chord  and  arc  lengths  are 
identical  for  the  range  of  arc  lengths  encountered. 

The  following  model  seems  appropriate  for  estimating  the 
error  in  distance.  Let  (X^,  Y^)  be  coordinates  of  the  signal  source  and 
(X^,  Y^)  be  the  coordinates  of  the  receiver.  Assume  further  that  the  errors 
in  placement  of  the  source  and  receiver  are  statistically  independent,  that 


(X^,  Yj)  and  (X^,  Y^)  are  samples  from  bivariate  normal  distributions  with 
identical  diagonal  covariance  matrices  (X  and  Y  errors  are  uncorrelated), 
and  that  the  mean  vectors  are  (o,  o),  and  (d,  o),  respectively.  In  general,  the 
distance  between  source  and  receiver  is  given  by 


d  =  V(x, .  x.r  +  (y,  -Y,r 


(A-  1 1) 


Since  X^,  X^,  Yj,  and  Y^  are  Gaussian  random  variables, 


it  is  seen  that 


X  =  (X2  -  Xt)  ~  N(d,  2o  ) 


(A-  12a) 


Y  =  (Y2  -  Y  t )  ~  N(o,  2o  ) 


(A-  12b) 


In  addition,  X  and  Y  are  uncorrelated.  This  latter  statement  follows  from 
the  diagonal-matrices  discussion  above.  The  distance,  which  is  not  normally 
distributed,  is  then 


D  =\/X2  +  Y2 


(A  -  1 3) 


It  should  be  recalled  that  the  assumptions  made  to  obtain  the 

A 

maximum- likelihood  estimate  b  for  the  slope  of  the  line  used  to  fit  the  data 

required  that  the  distance  error  be  normally  distributed.  It  appears  that  a 

2  2 

reasonable  approximation  can  be  made  for  the  distance  if  the  variance  o  <<  d  , 
since  under  this  condition  the  effects  of  errors  in  the  Y^  and  Y^  coordinates 
are  negligible  in  Equation  A- 13.  The  distance  is  modeled  as 


d  =  x2-x, 


(A  -  14) 


This  statement  can  be  written  as 


D  =  d  +  E 


(A-  15) 


A-4 


science  services  division 


where  E  is  a  Gaussian  random  variable  with  zero  mean  and  variance  of  2o  . 

x  2 
The  variance  in  source  or  receiver  placement  is  a  . 

Using  these  approximations,  the  development  of  the  estimate 

*  2 
b  is  valid,  provided  K  is  written  as 


(A-16) 


An  estimate  was  made  of  the  approximate  value  of  Equation  A-16 
using  previous  experience  as  a  guide.  The  value  obtained  from  this  procedure 
is  s  1.5  x  10  ^  sec^/km^  (a^  =  1.57  x  10  ^  sec^,  and  =  10.4  km^).  This 
value  is  used  in  all  final  computations  which  utilize  Equation  A-6  or  A-16. 

Once  the  line  segments  are  fitted  to  the  data  using  Equations 
A-2  and  A-8  on  each  segment,  standard  formulas^  are  employed  to  calculate 
the  value  of  each  parameter  of  the  model  using  the  fitted  lines,  slopes,  and 
intercept  values. 

At  this  point,  it  is  desirable  to  develop  an  expression,  or  at 
least  bounds,  for  the  variance  of  each  estimated  quantity.  The  ultimate  ob¬ 
jective  is  to  establish  confidence  limits  on  each  estimated  quantity.  An 

2 

approximate  expression  for  the  variance  a  (F)  of  a  function  F(X.,  X  ,  .  .  .  , 

17  t  2 

X  ) is  given  by 
n 

°2(r,  =  (l1)2a2<xi»+  (lr2)2  0 2  (x)2  +  •  •  • +  (if)2 "2 (x„'  ,A-17' 

when  the  jx^|  are  mutually  independent  and  includes  covariance  terms  of 
the  form 


(Si)  (S )  -  V 
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when  the  X.  are  not  mutually  independent.  The  quantities  for  which  estimates 
of  the  variance  are  desired  are  the  velocity  in  the  refractor,  the  dip  angle  of 
the  interface,  and  the  depths  to  the  refractor  at  each  end  of  the  reverse  pro¬ 
files.  The  velocity  in  the  refractor  is  given  by 


V.  = 


(A-  19) 


where  is  the  crustal  velocity,  assumed  to  be  without  error.  6  is  a  function 
of  the  estimates  of  the  slopes  of  the  fitted  lines  and  is  given  by 

9  *  j  |arc  sin  (v  ^  bdj  +  arc  sin  b^J  (A-20) 


where  b^  is  obtained  using  Equation  A-8  and  the  down-dip  traveltime  and  dis 
tance  data  and  where  b^  is  obtained  from  the  data  recorded  shooting  up-dip. 
The  dip  angle  between  the  crust  and  refractor  is 


r 

/  A  V 

A 

arc  sin 

(v.  bd) 

-  arc  sin 

(V.  bu)] 

(A  -  2  1) 


The  two  depths  are 


Vi  t 

1  ou 

A  A 

2  cos  (6)  cos  (5) 


(A-22) 


V  t  , 
1  od 


2  cos  (6)  cos  (6) 


(A-  23) 


where  t  and  t  ,  are  the  intercepts  of  the  line  fits  for  the  up-dip  and  down- 
ou  od 

dip  data,  respectively. 


In  the  succeeding  work,  approximate  values  for  the  variance 

of  each  of  these  quantities  are  obtained  assuming  the  covariance  terms  of 

Equation  A- 18  are  negligible  for  the  (t  ,  b.)  pairs.  This  follows  the  work  of 
15  01 

Davies.  Variances  of  the  estimated  quantities  related  to  the  fitting  of  the 
lines  to  the  data  sets  are^ 


2  A 
0  (b) 


(«'  • i!)' 

r  N  2  »  n  , 

£  (t  -  t)  -  b  £  (X.  -  X)  (t.  -  t) 

i=l  i=  1 

K2  (N  -  2) 

2  N  2  A  N 

K  £  (X.  -  X)‘  +  b  £  (X.  -  X)  (t.  - 

U1  i^l 

t) 

(A -24) 


a2  (.)  •  HiiL 

0  N 


+  (X)2  a2  <b)  +  b2  a2  (X) 


and 


(A-25) 


(A-26) 


(A-27) 


These  expressions  are  obtained  using  Equations  A- 17,  A-8, 
A- 10,  A-3,  and  A- 14.  It  can  now  be  shown  that  the  expressions  for  the 
desired  variance  are  given  by 


2  -  vi  cot  (9) 

a2  (V2>  =  — " — — — 

4  sin  (9) 


<bd> 

1  -  V^d 


a2  (bu) 

2  A2 

-  v: 

1  u 


2  *  V2  o2  (bd, 

a  (6)  = 


2  2  * 

V1  0  lbu' 


/,  2 

A  2\ 

4 

(>-vi 

b 

u/ 

(A-28) 


(A-29) 
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o  (h  \  = 
V  ou/ 


4Nu  cos2  (6)  cos2  (0) 


b  .  a2  (X  )  +  o2  ( Y  ) 
u  u  u' 


(A  -  30) 


Vi  1 

1  ou 


(bd  +  bu)\/'  -  vf  (b 


1  V  d/ 


tan  (S)  tan  (6)  +  tan2  (6)  >  a2  (b^ 


*  I  /  V  t  \ 

tan  (0)  —  /  1  ou 

bd * r  w-v,2 


tan  (6)  -  tan  (0) 


The  expression  for  g2^h  ^is  identical  to  Equation  A-30  with  the  subscript  u 
replaced  by  d  and  d  replaced  by  u  wherever  they  appear. 

The  previous  expressions,  obtained  using  Equation  A- 17,  are 
of  unknown  value  due  to  the  nonlinear  nature  of  Equations  A- 19  through  A-23. 
At  this  point,  sample  variance  will  be  considered  as  a  method  of  estimating 
the  true  population  variance. 

It  can  be  shown  that  the  sample  variance  is  an  optimum  esti¬ 
mate  in  the  sense  of  maximum  likelihood  of  the  true  variance  when  the 
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distribution  is  normal,  and  it  is  suspected  that  the  estimates  so  obtained 
will  be  good  estimates  provided  the  true  density  is  symmetrical  about  its 
peak  and  is  approximately  normal.  The  maximum-likelihood  estimator  just 
discussed  is  efficient  and  consistent,  and  it  possesses  a  property  called  in¬ 
variance.  For  these  reasons,  the  estimate  for  the  variance  of  each  quantity 
determined  by  Equations  A- 19,  A-21,  A-22,  and  A-23  will  be  the  correspond¬ 
ing  sample  variance.  Denotating  the  elements  |  V ,  6,  hQ  ,  h  by  {  m|  , 
the  unbiased  (normal  case)  estimate  of  the  variance. is 


— 5Z  (M.  -  M)2 
N-  1  i=1  » 


(A  -  3  1 ) 
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where  M,  the  sample  mean,  is  given  by 


N 

i=l 

However,  for  the  volume  of  data  available,  it  is  necessary  to  use  Equations 
A-28,  A-29,  and  A-30,  as  N  =  1  in  Equations  A-31  and  A-32. 

2 

Confidence  limits  can  now  be  determined  using  some  a  as  the 
variance  of  a  normal  distribution  of  zero  mean  and  confidence  level  N  such  as 


M. 


(A-32) 


A  =  M  +  P  (N)  o  (A-33a) 

and 

B  =  M  -  P  (N)  a  (A-33b) 


A  and  B  denote  the  upper  and  lower  N-percent  confidence  limits  on  the  com¬ 
puted  quantity  M  obtained  from  Equations  A- 19,  A-21,  A-22,  and  A-23.  P(N) 
is  the  standard-deviation  multiplier  obtained  for  N-percent  confidence  of  a 

normal  random  variable  of  zero  mean  and  unit  variance  (found  in  standard 

2 

tables);  a  is  obtained  from  o  .  The  computed  parameters  and  their  confi¬ 
dence  limits,  corr  puted  at  both  the  95-percent  |p(0.95)  =  1 . 96  J  and  the  90- 
percent  fp(0.90)  =  1 . 64 1  confidence  limits,  are  shown  in  Tables  A- 1  and  A-2. 


Since  the  a  priori  value  of  K  (Equation  A- 16)  indicated  that 
distance  errors  are  relatively  more  significant  than  time  errors,  a  one- 
variable-in-error  line  fit  for  each  set  of  data  points  was  also  considered. 

A 

For  X-in-error  line  fit,  the  slope  b^  is  given  by  Equations  A-8  and  A-9 
(when  K2  =  0)  as 


b 

x 


N 


(t  -  t)2 


i=  1 
N 


£  (X.  -  X)  (t.  -  t) 


(A-34) 


i=  1 
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The  variance  of  the  slope  estimate  is 


2 

a 


(A-35) 


where 


(X.  -  X)  (t.-t) 


( A-  36) 


The  average  value  of  a  estimated  using  Equation  2-36  was 

2  X  2 
approximately  16  km  ,  indicating  that  the  original  a  priori  value  of  a 

2  X 
(10.4  km  )  was  of  the  correct  order  of  magnitude  and  was  approximately 

correct. 


Essentially  the  same  values  occur  for  model  parameters  and 
the  model  parameters'  confidence  limits  when  it  is  assumed  that  both  time 
and  distance  are  subject  to  error  as  when  it  is  assumed  that  only  time  is 
subject  to  error  (Table  A-l).  It  may  be  concluded,  therefore,  that  no  sig¬ 
nificant  increase  in  accuracy  of  the  parameter  calculations  is  gained  by  using 
the  more  complicated  two-variables-in-error  line  fits  rather  than  standard 
least-square-error  line  fits.  It  also  may  be  concluded  that  distance-error 
effects  are  negligible  compared  to  timing-error  effects. 
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PARAMETERS  AND  CONFIDENCE  LIMITS,  SOUTHERN  REVERSE  PROFILE 
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